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Abstract 

Starting from Newton’s gravitational theory, we give a general in¬ 
troduction into the spherically symmetric solution of Einstein’s vac¬ 
uum field equation, the Schwarzschild(-Droste) solution, and into one 
specific stationary axially symmetric solution, the Kerr solution. The 
Schwarzschild solution is unique and its metric can be interpreted as 
the exterior gravitational field of a spherically symmetric mass. The 
Kerr solution is only unique if the multipole moments of its mass 
and its angular momentum take on prescribed values. Its metric can 
be interpreted as the exterior gravitational field of a suitably rotat¬ 
ing mass distribution. Both solutions describe objects exhibiting an 
event horizon, a frontier of no return. The corresponding notion of a 
black hole is explained to some extent. Eventually, we present some 
generalizations of the Kerr solution. 
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1 Prelude 1 

In Sec. 1.1, we provide some background material on Newton’s theory of grav¬ 
ity and, in Sec. 1.2, on the flat and gravity-free Minkowski space of special 
relativity theory. Both theories were superseded by Einstein’s gravitational 
theory, general relativity. In Sec. 1.3, we supply some machinery for formu¬ 
lating Einstein’s field equation without and with the cosmological constant. 

1.1 Newtonian gravity 

Newton’s gravitational theory is described—in particular tidal gravitational 
forces—and applied to a spherically symmetric body (a “star"). 


Gravity exists in all bodies universally and is proportional to the quan¬ 
tity of matter in each [■ ■ ■ ] If two globes gravitate towards each other, 
and their matter is homogeneous on all sides in regions that are equally 
distant from their centers, then the weight of either globe towards the 
other will be inversely as the square of the distance between the centers. 

Isaac Newton[136] (1687) 


The gravitational force of a point-like mass m 2 on a similar one of mass 
mi is given by Newton’s attraction law, 


2—>1 


= -G 


m 1 m 2 r 


( 1 ) 


where G is Newton’s gravitational constant (CODATA 2010), 


G = 6.67384(80) x lfT 11 . 

The vector r := ri — r 2 points from m 2 to mi, see the Fig. 1. 

According to actio = reactio (Newton’s 3rd law), we have F 2 _>.i = — Fi_> 2 . 
Thus a complete symmetry exists of the gravitational interaction of the 
two masses onto each other. Let us now distinguish the mass m 2 as field¬ 
generating active gravitational mass and mi as (point-like) passive test- 
mass. Accordingly, we introduce a hypothetical gravitational field as describ¬ 
ing the force per unit mass (m 2 M, m ): 


F _ GM r 
m Ir1 2 |r 


( 2 ) 


With this definition, the force acting on the test-mass m is equal to field 
Warts of Secs.l & 2 are adapted from our presentation[83] in Falcke et al.[56]. 
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Figure 1: Two mass points m\ and m 2 attracting each other in 3-dimensional 
space, Cartesian coordinates x , y, z. 



Figure 2: The “source” M attracts the test mass m. 
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strength x gravitational charge (mass) or F^ m — mi, in analogy to elec¬ 
trodynamics. The active gravitational mass M is thought to emanate a 
gravitational field which is always directed to the center of M and has the 
same magnitude on every sphere with M as center, see Fig. 2. Let us now in¬ 
vestigate the properties of the gravitational field (2). Obviously, there exists 
a potential 

M 

0 = -G —, f = -V0. (3) 

|r| 

Accordingly, the gravitational field is curl-free: V X f = 0. 

By assumption it is clear that the source of the gravitational field is the 
mass M. We find, indeed, 


V-f = -47rGMh 3 (r), (4) 

where h 3 (r) is the 3-dimensional (3d) delta function. By means of the Laplace 
operator A := V-V, we infer for the gravitational potential 

A 0 = 47t G M S 3 (r). (5) 

The term M h 3 (r) may be viewed as the mass density of a point mass. Eq.(5) 
is a 2nd order linear partial differential equation for 0. Thus the gravitational 
potential generated by several point masses is simply the linear superposition 
of the respective single potentials. Hence we can generalize the Poisson 
equation (5) straightforwardly to a continuous matter distribution p( r): 

A 0 = 47 t G p. (6) 

This equation interrelates the source p of the gravitational field with the 
gravitational potential 0 and thus completes the quasi-field theoretical de¬ 
scription of Newton’s gravitational theory. 

We speak here of quasi -field theoretical because the field 0 as such rep¬ 
resents a convenient concept. However, it has no dynamical properties, no 
genuine degrees of freedom. The Newtonian gravitational theory is an action 
at a distance theory (also called mass-interaction theory ). When we remove 
the source, the field vanishes instantaneously. Newton himself was very un¬ 
happy about this consequence. Therefore, he emphasized the preliminary 
and purely descriptive character of his theory. But before we liberate the 
gravitational field from this constraint by equipping it with its own degrees 
of freedom within the framework of general relativity theory, we turn to some 
properties of the Newtonian theory. 

A very peculiar fact characteristic to the gravitational field is that the 
acceleration of a freely falling test-body does not depend on the mass of 
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Figure 3: Tidal forces emerging between two freely falling particles and de¬ 
forming a spherical body. 


this body but only on its position within the gravitational held. This comes 
about because of the equality (in suitable units) of the gravitational and the 
inertial mass: 

inertial .. grav _ 

m r = F = m f . 

This equality has been well tested since Galileo’s time by means of pendulum 
and other experiments with an ever increasing accuracy, see Will [189]. 

In order to allow for a more detailed description of the structure of a 
gravitational held, we introduce the concept of tidal force. This can be best 
illustrated by means of Fig. 3. In a spherically symmetric gravitational held, 
for example, two test-masses will fall radially towards the center and thereby 
get closer and closer. Similarly, a spherical drop of water is deformed to an 
ellipsoidal shape because the gravitational force at its bottom is bigger than 
at its top, which has a greater distance to the source. If the distance between 
two freely falling test masses is relatively small, we can derive an explicit 
expression for their relative acceleration by means of a Taylor expansion. 
Consider two mass points with position vectors r and r + hr , with |hr| <C 1. 
Then the relative acceleration reads 


ha = [f(r + hr) — f(r)] = hr • (V/). (8) 

We may rewrite this according to (the sign is conventional, d/dx a =: d a , 
x 1= x, x 2 = y , x 3 = z ) 

K ab :=-(Vf) ab = -d a f b , a, 5=1,2,3. (9) 
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We call K ab the tidal force matrix. The vanishing curl of the gravitational 
held is equivalent to its symmetry, K ab = K ba . Furthermore, K ab = d a d b (p. 
Thus, the Poisson equation becomes, 

3 

K aa = trace K = Air G p. (10) 

a =1 

Accordingly, in vacuum K ab is trace-free. 

Let us now investigate the gravitational potential of a homogeneous star 
with constant mass density p Q and total mass M & = (4/3) 7r p 0 . For 
our Sun, the radius is R@ = 6.9598 x 10 8 m and the total mass is M = 
1.989 x 10 30 kg. 

Outside the sun (in the idealized picture we are using here), we have 
vacuum. Accordingly, p( r) = 0 for |r| > R@. Then the Poisson equation 
reduces to the Laplace equation 

A0 = O, for r > R q . (11) 

In 3d polar coordinates, the r-dependent part of the Laplacian has the form 
(1/r 2 ) d r (r 2 d r ). Thus (11) has the solution 

4>=~ + P, ( 12 ) 

r 

where a and /3 are integration constants. Requiring that the potential tends 
to zero as r goes to infinity, we get f3 = 0. The integration constant a will 
be determined from the requirement that the force should change smoothly 
as we cross the star’s surface, that is, the interior and exterior potential and 
their first derivatives have to be matched continuously at r = R & . 

Inside the star we have to solve 


A f> = 4:7: G p Q , for r < R Q . (13) 

We find 

4> — q nGpQ r 2 H-h C *2 , (14) 

3 r 

with integration constants C\ and C^- We demand that the potential in the 
center r = 0 has a finite value, say <f>Q. This requires Ci—0. Thus 


, 2 9 ^ GM(r ) , 

0 =-7rGp Q r-+ 0 O = —2r +< ^ 0 ’ 

where we introduced the mass function M{r ) = (4/3) tt r 3 p e which measures 
the total mass inside a sphere of radius r. 
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Figure 4: Newtonian potential of a homogeneous star. 


Continuous matching of 0 and its first derivatives at r = finally yields: 

_ ^ "Trf f° r l r l — ’ 

= ^ 1\/T ‘in A/T 


G 


M ra , ,2 3 GMp 

2R% |r| 2 R q 


© 


© 


The slope of this curve indicates the magnitude of the gravitational force, the 
curvature (2nd derivative) the magnitude of the tidal force (or acceleration). 








1.2 Minkowski space 

When, in a physical experiment, gravity can be safely neglected, we seem to 
live in the flat Minkowski space of special relativity theory. We introduce 
the metric of the Minkowksi space and rewrite it in terms of so-called null 
coordinates, that is, we use light rays for a parametrization of Minkowski 
space. 


Henceforth space by itself, and time by itself, are doomed to 
fade away into mere shadows, and only a kind of union of 
the two will preserve an independent reality. 

Hermann Minkowski (1908) 

It was Minkowski who welded space and time together into spacetime, 
thereby abandoning the observer-independent meaning of spatial and tem¬ 
poral distances. Instead, the spatio-temporal distance, the line element, 

ds 2 = —c 2 d.t 2 + dx 2 + dy 2 + dz 2 

is distinguished as the invariant measure of spacetime. The Poincare (or 
inhomogeneous Lorentz) transformations form the invariance group of this 
spacetime metric. The principle of the constancy of the speed of light is 
embodied in the equation ds 2 = 0. Suppressing one spatial dimension, the 
solutions of this equation can be regarded as a double cone. This light cone 
visualizes the paths of all possible light rays arriving at or emitted from 
the cone’s apex. Picturing the light cone structure, and thereby the causal 
properties of spacetime, will be our method for analyzing the meaning of the 
Schwarzschild and the Kerr solution. 

Null coordinates 

We first introduce so-called null coordinates. The Minkowski metric (with 
c = 1), in spherical polar coordinates reads 

ds 2 = — dt 2 + dr 2 + r 2 (dd 2 + sin 2 6 dfl 2 ) = —d.t 2 + dr 2 + r 2 dM 2 . (17) 

We define advanced and retarded null coordinates according to 

v := t + r , u t — r , (18) 

and find 

ds 2 = —dv du + - (v — u) 2 dM 2 . (19) 
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Figure 5: Minkowski spacetime in null coordinates 


In Fig. 5 we show the Minkowski spacetime in terms of the new coordinates. 
Incoming photons, that is, point-like particles with velocity r = —c = — 1, 
move on paths with v = const. Correspondingly, we have for outgoing pho¬ 
tons u = const. The special relativistic wave-equation is solved by any func¬ 
tion f(u) and f{v). The surfaces f{u) = const, and f{v ) = const, represent 
the wavefronts which evolve with the velocity of light. The trajectory of ev¬ 
ery material particle with r < c — 1 has to remain inside the region defined 
by the surface r = t. In an (r, t)-diagram this surface is represented by a 
cone, the so-called light cone. Any point in the future light cone r = t can be 
reached by a particle or signal with a velocity less than c. A given spacetime 
point P can be reached by a particle or signal from the spacetime region 
enclosed by the past light cone r = —t. 
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Penrose diagram 

We can map, following Penrose, the infinitely distant points of spacetime into 
finite regions by means of a conformal transformation which leaves the light 
cones intact. Then we can display the whole infinite Minkowski spacetime 
on a (finite) piece of paper. Accordingly, introduce the new coordinates 

v := arctan v , u := arctan u , for — 7r/2 < (v , u ) < +7t/2 . (20) 


Then the metric reads 


ds 2 = 


cos 2 v cos 2 u 


—dv du -\— sin 2 (v — u ) dVt 2 
4 v ' 


( 21 ) 


We can go back to time- and space-like coordinates by means of the trans¬ 
formation 

t v + u, f := v — u , (22) 

see (18). Then the metric reads, 

, 2 — d,t 2 + dr 2 + sin 2 f dQ 2 

ds 2 = -=-=- 

4 cos 2 ^ cos 2 ^ 

that is, up to the function in the denominator, it appears as a flat metric. 
Such a metric is called conformally flat (it is conformal to a static Einstein 
cosmos). The back-transformation to our good old Minkowski coordinates 
reads 



1 

f t + r 


t — r\ 

(24) 

f “ 2 

( tan 2 

+ tan 

2 )• 

1 

( t + f 


t — r\ 

(25) 

r — - 

2 

(tan 2 

— tan 

2 ■ 


Our new coordinates t, f extend only over a finite range of values, as can 
be seen from (24) and (25). Thus, in the Penrose diagram of a Minkowski 
spacetime, see Fig. 6, we can depict the whole Minkowski spacetime, with 
a coordinate singularity along r = 0. All trajectories of uniformly moving 
particles (with velocity smaller than c) emerge form one single point, past 
infinity I~, and all will eventually arrive at the one single point J + , namely 
at future infinity. All incoming photons have their origin on the segment 
Z~ (script I~ or “scri minus”), light-like past-infinity, and will run into the 
coordinate singularity on the t-axis . All outgoing photons arise from the 
coordinate singularity and cease on the line Z + , light-like future infinity (“scri 
plus”). The entire spacelike infinity is mapped into the single point 7°. For 
later reference we collect these notions in a table: 


11 











t = const 



r 


Figure 6: Penrose diagram of Minkowski spacetime. 


Table 1. The different infinities in Penrose diagrams 

/- 

timelike past inhnity 

origin of all particles 

1+ 

timelike future infinity 

destination of all particles 

h 

spacelike infinity 

inaccessible for all particles 

x~ 

lightlike past infinity 

origin of all light rays 

x + 

lightlike future infinity 

destination of all light rays 


Now, we have a really compact picture of the the Minkowski space. Next, 
we would like to proceed along similar lines in order to obtain an analogous 
picture for the Schwarzschild spacetime. 
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1.3 Einstein’s field equation 

We display our notations and conventions for the differential geometric tools 
used to formulate Einstein’s field equation. 


We assume that our readers know at least the rudiments of general relativ¬ 
ity (GR) as represented, for instance, in Einstein’s Meaning of Relativity[ 50], 
which we still recommend as a gentle introduction into GR. More advanced 
readers may then want to turn to Rindler[165] and/or to Landau-Lifshitz[106]. 

We assume a 4d Riemannian spacetime with (Minkowski-)Lorentz signa¬ 
ture (-b ++), see Misner, Thorne, and Wheeler[ ]. Thus, the metric 

field, in arbitrary holonomic coordinates x* 1 , with p = 0,1, 2, 3, reads 

g = ds 2 = g jW dx M ® dx u . (26) 


By partial differentiation of the metric, we can calculate the Christoffel sym¬ 
bols (Levi-Civita connection) 

:= ^ g^ {d a g Pl + d p g ia - d^g a y) . (27) 


This empowers us to determine the geodesics (curves of extremal length) of 
the Riemannian spacetime: 


D 2 x a d 2 x a dx 11 dx v 

Dt 2 dr 2 l ‘ v dr dr 


(28) 


This equation can be read as a vanishing of the 4d covariant acceleration. If 
we define the 4-velocity u a := dx a /dr , then the geodesics can be rewritten 
as 


Du c 


du c 


+ T“ u y/ = 0 . 


Dr dr ' ~ (29) 

In a neighborhood of any given point in spacetime we can introduce 
Riemannian normal coordinates, which are such that the Christoffels van¬ 
ish at that point. In order to find a tensorial measure of the gravitational 
held, we have to go one differentiation order higher. By partial differentiation 
of the Christoffels, we find the Riemann curvature tensor 2 


R^vap ■— 2 {d [a T% m + W a[a T a \um) ■ (30) 

2 Always symmetrizing of indices is denoted by parentheses, (a/3) := {a/3 + /3a}/ 2!, 
antisymmetrization by brackets [a/3] := { a d ~ a /3}/2!, with corresponding generalizations 
{a/3y) := {+a/3y+(3ya+'ya/3 + - ■ ■ }/3!, etc.; indices standing between two vertical strokes 
| are excluded from the (anti)symmetrization process, see Schoutenf 70]. 
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The curvature is doubly antisymmetric, its two index pairs commute, and its 
totally antisymmetric piece vanishes: 

0 j Rfiu(aP) 0 ■ R/iuap Rapfiv ■ R[fiuaP] 0 ■ (31) 

If we define collective indices A,B,.. = 1,..., 6 for the antisymmetric index 
pairs according to the rule {01,02,03; 23,31,12} —> {1,2, 3; 4, 5, 6 }, then 
the algebraic symmetries of (31) can be rephrased as 

Rab = Rba , trace(i? J 4e) = 0. (32) 

Thus, in 4d the curvature can be represented as a trace-free symmetric 6 x 6 - 
rnatrix. Hence it has 20 independent components. 

With the curvature tensor, we found a tensorial measure for the grav¬ 
itational field. Freely falling particles move along geodesics of Riemannian 
spacetime. What about the tidal accelerations between two freely falling par¬ 
ticles? Let the “infinitesimal” vector n a describe the distance between two 
particles moving on adjacent geodesics. A standard calculation! 7], linear 
to the order of n, yields the geodesic deviation equation 

(33) 

This equation describes the relative acceleration of neighboring particles, 
similar as ( 8 ) and (9) in the Newtonian case. The role of the tidal matrix 
K ab is taken over by K. a s ■— u^u 1 R a 

By contraction of the curvature, we can define the 2nd rank Ricci tensor 
R^v and the curvature scalar R, respectively: 

RnV := R a na,V , R := g^R^iu • (34) 

For convenience, we can also introduce the Einstein tensor G^ := R^ — 
\QhvR- The curvature with its 20 independent components can be irreducibly 
decomposed into smaller pieces according to 20 = 10 + 9 + 1. The Weyl 
curvature tensor C a p lb is trace-free and has 10 independent components, 
whereas the trace-free Ricci tensor has 9 components and the curvature scalar 
just 1 . 

Now we have all the tools for displaying Einstein’s field equation. With G 
as Newton’s gravitational constant and c as velocity of light, we define Ein¬ 
stein’s gravitational constant k := 87 rG/c 4 . Then, the Einstein field equation 
with cosmological constant A reads 

R^lv 2 9/ivR t ^g»v ^ ■ (35) 
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The source on the right-hand side is the energy-momentum tensor of matter. 
The vacuum field equation, without cosmological constant, simply reduces to 
R /1L , = 0. Mostly this equation will keep us busy in this article. A vanishing 
Ricci tensor implies that only the Weyl curvature C a ^s 7^ 0. Accordingly, 
the vacuum field in GR (without A) is represented by the Weyl tensor. 

Eq.(35) represents a generalization of the Poisson equation (10). There, 
the contraction of the tidal matrix is proportional to the mass density; in 
GR, the contraction of the curvature tensor is proportional to the energy- 
momentum tensor. 

The physical mass is denoted by M. Usually, we use the mass parameter, 
m := (= pf ■ The Schwarzschild radius reads rs '■= 2 m = =^r~. Usually we put 
c = 1 and G — 1. We make explicitly use of G and c as soon as we stress 
analogies to Newtonian gravity or allude to observational data. 
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2 The Schwarzschild metric (1916) 

Spatial spherical symmetry is assumed and a corresponding exact solution 
for Einsteins theory searched for. After a historical outline (Sec.2.1), we 
apply the equivalence principle to a freely falling particle and try to imple¬ 
ment that on top of the Minkowskian line element. In this way, we heuristi- 
cally arrive at the Schwarzschild metric (Sec.2.2). In Sec.2.3, we display the 
Schwarzschild metric in six different classical coordinate systems. We out¬ 
line the concept of a Schwarzschild black hole in Sec. 2. f. In Secs.2.5 and 2.6, 
we construct the Penrose diagram for the Schwarzs child (-Kruskal) spacetime. 
We add electric charge to the Schwarzschild solution in Sec.2.7. The interior 
Schwarzschild metric, with matter, is addressed in Sec.2.8. 


It is quite a wonderful thing that from such an abstract 
idea the explanation of the Mercury anomaly emerges so in¬ 
evitably. 

Karl Schwarzschild[171] (1915) 


2.1 Historical remarks 

The genesis of the Schwarzschild solution (1915/16) is described. In par¬ 
ticular, we show that Droste, a bit later than Schwarzschild, arrived at the 
Schwarzschild metric independently. He put the Schwarzschild solution into 
that form in which we use it today. 

The first exact solution of Einstein’s field equation was born in hospi¬ 
tal. Unfortunately, the circumstances were more tragic than joyful. The 
astronomer Karl Schwarzschild joined the German army right at the begin¬ 
ning of World War I and served in Belgium, France, and Russia. At the 
end of the year 1915, he was admitted to hospital with an acute skin dis¬ 
ease. There, not far from the Russian front, enduring the distant gunfire, 
he found time to “stroll through the land of ideas” of Einstein’s theory, 
as he puts it in a letter to Einstein 3 dated 22 December 1915. According 
to this letter, Schwarzschild started out from the approximate solution in 
Einstein’s “perihelion paper”, published November 25th. Since presumably 
letters from Berlin to the Russian front took a few days, Schwarzschild [17: 
found the solution within about a fortnight. Fortunately, the premature field 
equation of the “perihelion paper” is correct in the vacuum case treated by 
Schwarzschild. 

3 The letters from and to Einstein can be found in Einstein’s Collected Works[51], see 
also Schwarzschild’s Collected Works[171]. 
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In February 1916, Schwarzschild[173] submitted the spherically symmet¬ 
ric solution with matter—the “interior Schwarzschild solution”—now based 
on Einstein’s final field equation. In March 1916, he was sent home were he 
passed away on 11 May 1916. 

The field equation used by Schwarzschild requires detg = — 1. To fulfill 
this condition, he uses modified polar coordinates (Schwarzschild’s original 
notation used), 

r 3 

Xi — — , x 2 = — cos 9 , x 3 = 0 , X4 — t. 

The spherically symmetric ansatz then reads 

dr 2 

ds 2 = U dx 4 - /1 dx{ - f 2 ' % - f 3 dxl (1 - xf) , 

1 -*2 

where /1 to A are functions of X\ only. The solution turns out to be 

A = ^4 1 _^ a / R > A = A = R 2 , / 4 = 1 - a/-R, i? = (r 3 + a 3 ) 1/3 . 

In this article, as well as in his letter to Einstein, he eventually returns to 
the usual spherical polar coordinates, 

ds 2 = (1 — a/R) dt 2 - ( ^' — R 2 ( dO 2 + sin 2 6 def) 2 ) , R = (r 3 + a 3 ) 1//3 . 

1 — a/R 

This looks like the Schwarzschild metric we are familiar with. One should 
note, however, that the singularity at R = a is (as we know today) a coordi¬ 
nate singularity, it corresponds to r = 0. In the early discussion the meaning 
of such a singularity was rather obscure. Flamm[60] in his 1916 article on 
embedding constant time slices of the Schwarzschild metric into Euclidean 
space mentions “the oddity that a point mass has an finite circumference of 
27ra”. 

In 1917, Weyl[188] talks of the “inside” and “outside” of the point mass 
and states that “in nature, evidently, only that piece of the solution is real¬ 
ized which does not touch the singular sphere.” In Hilbert’s[ 86 ] opinion, the 
singularity R = a indicates the illusiveness of the concept of a pointlike mass; 
a point mass is just the limiting case of a spherically symmetric mass dis¬ 
tribution. Illuminating the interior of “Schwarzschild’s sphere” took quite a 
while and it was the discovery of new coordinates which brought first elucida¬ 
tions. Lanczos[104], in 1922, clearly speaks out that singularities of the met¬ 
ric components do not necessarily have physical significance since they may 
vanish in appropriate coordinates. However, it took another 38 years to find 
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a maximally extended fully regular coordinate system for the Schwarzschild 
metric. We will become acquainted with these Kruskal/Szekeres coordinates 
in Sec. 2.5. 

Schwarzschild’s solution, published in the widely read minutes of the 
Prussian Academy, communicated by Einstein himself, nearly instantly trig¬ 
gered further investigations of the gravitational held of a point mass. Already 
in March 1916, Reissner[ .63], a civil engineer by education, published a gen¬ 
eralization of the Schwarzschild metric, including an electrical charge; this 
was later completed by Weyl[188] and by Nordstrom[ 10]. Today it is called 
Reissner-Nordstrom solution. 

Nevertheless, one should not ignore the Dutch twin of Schwarzschild’s 
solution. On 27 May 1916, Drostef ] communicated his results on “the 
held of a single centre in Einstein’s theory of gravitation, and the motion 
of a particle in that held” to the Dutch Academy of Sciences. He presents 
a very clear and easy to read derivation of the metric and gives a quite 
comprehensive analysis of the motion of a point particle. Since 1913, he had 
been working on general relativity under the supervision of Lorentz at Leiden 
University. Published in Dutch, Droste’s results are fairly unknown today. 
Einstein, probably informed by his close friend Ehrenfest, rather appreciated 
Droste’s work, praising the graceful mathematical style. Weyl[ 188] also cites 
Droste, but in Hilbert’s[86] second communication the reference is not found. 
Einstein, Hilbert, and Weyl always allude to “Schwarzschild’s solution”. 

After Droste took his PhD in 1916, he worked as school teacher and 
eventually became professor for mathematics in Leiden. He never resumed his 
work on Einstein’s theory and his name faded from the relativistic memoirs. 
In Leiden, people like Lorentz, de Sitter, Nordstrom, or Fokker learned about 
the gravitational held of a point mass primarily from Droste’s work. Thus, 
the name “Schwarzschild-Droste solution” would be quite justified from a 
historical point of view. 

The importance of the Schwarzschild metric is made evident by the Birk- 
hoff[l6] theorem 4 : For vanishing cosmological constant, the unique spheri¬ 
cally symmetric vacuum spacetime is the Schwarzschild solution, which can 
be expressed most conveniently in Schwarzschild coordinates, see Table 3, 
entry 1. Thus, a spherically symmetric body is static (outside the horizon). 
In particular, it cannot emit gravitational radiation. Moreover, the asymp¬ 
totic Minkowskian behavior of the Schwarzschild solution is dictated by the 
solution itself, it is not imposed from the outside. 

4 The “Birkhoff” theorem was discovered by Jebsen[95], BirkhofffUi], and 
Alexandrow[4]. For more details on Jebsen, see Johansen & Ravndal[96]. The objections 
of Ehlers & Krasinski[48] appear to us as nitpicking. 
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2.2 Approaching the Schwarzschild metric 

We start from an ansatz for the metric of an accelerated motion in the radial 
direction and combine it, in the sense of the equivalence principle, with the 
free-fall velocity of a particle in a Newtonian gravitational field. In this way, 
we find a curved metric that, after a coordinate transformation, turns out to 
be the Schwarzschild metric. 

Einstein, in his 1907 Jahrbuch article [49], suggests the generalization of 
the relativity principle to arbitrarily accelerated reference frames. 

A plausible notion of a (local) rest frame in general relativity is a frame 
where the coordinate time is equal to the proper time (for an observer spa¬ 
tially at rest, of course). For a purely radial motion, the following metric 
would be an obvious ansatz, see also Visser[185]: 

ds 2 = — dt 2 + [dr + /(r) dt] 2 + r 2 dfl 2 , with dfl 2 := dd 2 + sin 2 9d(j) 2 . 

For d<f> = 0, d6 = 0, and dr/dt = —f(r), we have ds 2 = —dt 2 . Thereby, — f(r ) 
is identified as a kind of “radial infall velocity”. Note also that constant time- 
slices, dt = 0, are Euclidean. 

In Newtonian gravity, a particle falling from infinity towards the origin 
picks up a velocity 


dr / / 2 GM 

s =« = -^*w = -d— 


1 2 / \ if,/ \ GM 

-mv (r) = m<P(r) = m —— 


(36) 

Here, $ is the absolute value of the Newtonian potential of a spherical body 
with mass M. 

Hence, in some Newtonian limit, we demand f(r) —» \/2$. This leads to 
the metric 


ds 2 = — dt 2 + (dr + y/2^> dt ) 2 + r 2 dQ 2 , (37) 


where we allow for an arbitrary potential if = if{r). This metric generates 
curvature. The calculations can be conveniently done even by hand. The 
Ricci tensor reads 


R 0 ° = Ri 1 = - d r d r (r fj) = 0 , 
r 


TD 2 _ o 3 _ 

/to — ill — 


2 d r (rif) 


= 0 


The equations i?o° = 0 = Ri 1 are mere integrability conditions of the R 2 2 = 
0 = R 3 3 relations. Hence, rif is determined by its first order approximation 
alone and reads 


i> = 


a 

1 

r 
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with a as an unknown constant so far. By construction, we have 


dr 

dt 


-Vw = 



2 GM 


a = GM —: m,. 


The metric (37), expanding the parenthesis and collecting the terms in 
front of dt 2 , reads 


, , ( 2 GM\ , 9 I2GM , , , 9 9 ~ 

ds 2 — — I 1- J dt 2 + 2 J —-- dtdr + dr 2 + r 2 dfl 2 . (38) 

Using different methods, this metric was derived by Gullstrand[70] in May 
1921. Gullstrand claimed to have found a new spherically symmetric solution 
of Einstein’s field equation. In his opinion 5 , this showed the ambiguity of 
Einstein’s field equation. However, the metric is of the form 

,9 , , 9 , , , 9 9 „ 2 GM hGM 

ds 2 = —Adt 2 + 2Bdtdr + dr 2 + r 2 dfl 2 , A:— 1-, B := \ -, 

r V r 

and can be diagonalized by completing the square via 

— ^ dr^j + ^1 + d r2 + r2 ■ 


ds 2 = —A ( dt 


Introducing a new time coordinate, 


dt s 


dt — 


— dr 
A 


or, explicitly, t s = t - (zr^^r- - 4 GM Artanh, we arrive at (A 
and B re-substituted) 


ds 2 



dt 2 s + 


/ _ 2 GM 

y r 


-l 

dr 2 + r 2 dQ 2 . 


In contrast to what Gullstrand was aiming at, he “just” rederived the Schwarz- 
schild metric. 

Later, applying a coordinate transformation to the Schwarzschild met¬ 
ric, Painleve[ ] obtained the metric (38) independently and presented his 
result in October 1921. His aim was to demonstrate the vacuity of ds 2 by 

5 Gullstrand, who was a member of the Nobel committee, was responsible that Einstein 
did not get his Nobel prize for relativity theory. He thought that GR is untenable. 
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showing that an exact solution does not determine the physical geometry 
and is therefore meaningless. In a letter (Dec. 7 th 1921) to Painleve, Ein¬ 
stein stresses on the contrary the meaninglessness of the coordinates! In the 
words of Einstein himself (our translation): “.. .merely results obtained by 
eliminating the coordinate dependence can claim an objective meaning.” 

In the subsequent section, we will meet the Schwarzschild metric in many 
different coordinate systems. All of them have their merits and their short¬ 
comings. 

Using Gullstrand-Painleve coordinates for the Schwarzschild metric does 
not change the physics, of course. However, as a coordinate system it is what 
Gustav Mie[126] calls a sensible coordinate system. In contrast to many 
other coordinate systems, the physics looks quite like we are used to. As an 
example, we analyze the motion of a radial infalling particle in Schwarzschild 
and Gullstrand-Painleve coordinates. 

The equations of motion for point particles in general relativity are ob¬ 
tained via the geodesic equation (28). It can be shown that this equa¬ 
tion is equivalent to the solution of the variational principle 5 f a ds 2 = 
5 f x a x 13 g a p dr 2 . We choose the proper time r for the parametrization of 
the curve, the dot denotes the derivative with respect to r. In the present con¬ 
text, we are only interested in the velocity of particles along ingoing geodesics 
(“freely falling particles”). For time-like geodesics we have — 1 = This 
allows the algebraic determination of r provided we know t. Since we con¬ 
sider static metrics here, t is a cyclic variable and ^ = K = const. 

The constant is determined by the boundary condition r = 0 for r —y oo. 
The calculation yields: 


Table 2. Velocities in different coordinate systems 6 




Schwarzschild 

Gullstrand-Painleve 


coordinate 

_|_ ^ 2GM ^ / 2GM 



velocity 

2GM 

V r 

particles 

dt 

proper 




velocity 

dr 

dr 




coordinate 



light rays 

velocity 

dr 

d,t 

± (1_ 2^M) 

±i AV 


6 The velocities of outgoing particles are valid only for the boundary condition specified. 
The coordinate velocity for outgoing particles in GP coordinates does not fit in our table 
and is thus suppressed. 
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The difference between the coordinate systems appears in the first line 
of Table 2: In Gullstrand-Painleve coordinates, the coordinate velocity of a 
freely infalling particle increases smoothly towards the center. Nothing spe¬ 
cial happens at r = 2 GM. From a given position, the particle will plunge into 
the center in a finite time. Even numerically this looks quite Newtonian. In 
contrast, the velocity with respect to Schwarzschild coordinates approaches 
zero as the particle approaches r = 2 GM. Hence, the particle apparently 
will not be able to go further than r = 2 GM. 

For the Gullstrand-Painleve metric for incoming light the radial coordi¬ 
nate velocity is always larger in magnitude than —1, at r = 2 GM it is —2, 
for outgoing rays it vanishes at r = 2 GM and is negative for r < 2 GM. 

Taking the mere numerical values is misleading. Contemplate for incom¬ 
ing light 

(—) 

\ dt / particle _ 

—) ~ l + 

dt) light 

So the particle is always slower than light, however it approaches the velocity 
of light when approaching r = 0. 

The Gullstrand-Painleve form of the metric is regular at the surface r = 
2 GM. This shows that it is not any kind of barrier, but this observation was 
not made until much later, see Eisenstaedt [52]. 

2.3 Six classical representations of the Schwarzschild 
metric 

As we mentioned, a coordinate system should be chosen according to its con¬ 
venience for describing a certain situation. In the following table (Table 3), 
we collect six widely used forms of the Schwarzschild metric. 




< i 


2GM 
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2.4 The concept of a Schwarzschild black hole 


We first draw a simple picture of a black hole. The event horizon and the 
stationary limit emerge as characteristic features. These are subsequently 
defined in a more mathematical way. 

In 1783 John Michell communicated his thoughts on the means of dis¬ 
covering the Distance, magnitude, etc. of the fixed stars, in consequence of 
the diminuation of the velocity of their light ... [125] to the Royal Society 
in London. In the context of Newton’s particle theory of light, he calcu¬ 
lated that sufficiently massive stars exhibit a gravitational attraction to such 
vast an amount that even light could not escape. A few years later (1796) 
Pierre-Simon Laplace published similar ideas. 

In modern notation, we may reconstruct the arguments as follows. We 
throw a mass m from the surface of the Earth, assuming that there were no 
air, in upward direction with an initial velocity v. It will always fall back, 
unless its initial velocity reaches a sufficiently high value ry sca pe providing 
the mass with such a kinetic energy that it can overpower the gravitational 
attraction of the Earth. Energy conservation yields then immediately the 
formula 

h GM @ 

^escape \ / 1 

y -Re 

where G is Newton’s gravitational constant and M and the mass and the 
radius of the spherically conceived Earth, respectively. 

For the Earth we find ^escape ~ 11.2 km/s. If we now compress the Earth 
appreciably (thought experiment!) until the escape velocity coincides with 
the speed of light u esC ape = c, its compressed “Schwarzschild” radius becomes 
r e = 2 GM^/c 2 ~ 1 cm. For the Sun, with its mass M 0 , we have' 


r& 


2 GM q 


pc 3 km. 


At any smaller radius the light will be confined to the corresponding body. 
This is an intuitive picture of a spherically symmetric invisible “black hole” . * 8 

It is very intriguing to see how far-sighted Michell anticipated the status 
of today’s observational black hole physics: 

' For the sake of clarity, we display here the speed of light c explicitly. 

8 The phrase “black hole” is commonly associated with Wheeler (1968). It appears 
definitely earlier in the literature: In the January 1964 edition of the Science News Letter 
the journalist Ann Ewing entitled her report at the meeting of the American Association 
for the Advancement of Science in Cleveland “Black Holes” in space. And if you have 
a look into an arbitrary English language dictionary published before ca. 1970, you will 
learn that “black hole” refers to a notorious dungeon in Calcutta (now Kolkata) in the 
18th century, apparently a place of no return ... 
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If there should really exist in nature bodies, whose density is not less than that 
of the sun, and whose diameters are more than 500 times the diameter of the sun, 
since their light could not arrive at us; [... ] we could have no information from 
sight; yet, if any other luminous bodies should happen to revolve about them we 
might still perhaps from the motions of these revolving bodies infer the existence 
of the central ones with some degree of probability ... 

This could be a verdict on the current observations of the black hole Sgr 
A* (“Sagittarius A-star”) in the center of our Milky Way—and this is not a 
thought experiment—for a popular account, see Sanders[1 ]. Sgr A* has a 

mass of about 4 x 10 6 M 0 . Thus its Schwarzschild radius is far from being 
minute, it is about 3 x 4 x 10 6 km or about 17 solar radii. 



Figure 7: Not quite seriously: “Schwarzschild” (left) versus “Kerr” (right) 

A cautionary remark has to be made, though, see Penrose[150]. In New¬ 
tonian gravity c has no absolute meaning like in special relativity. It is con¬ 
ceivable that the speed of light in strong Newtonian gravitational fields could 
be larger than c. Consequently, the Michcll type argument becomes only per¬ 
tinent if c is the maximal speed for all phenomena like in the Minkowski space 
of special relativity, or, if gravity is involved, in the Riemannian space of GR. 

Let us follow the way of visualizing the black hole concept by means of 
everyday physics a bit further: We explore the Schwarzschild and, later in 
Sec. 3.4, the Kerr spacetime by boat. Schwarzschild spacetime is mimicked 
by a hole in a lake in which the surrounding water plunges simply radially 
without whirling around (Monticello Dam, California). The water flowing 
towards the hole will drag our boat to the center. Our boat may move 
around quite freely as long as the current is weak. 

However, at some distance from the hole, the current becomes so strong 
that our boat, engines working at their maximum power, merely can keep its 
position. This is the stationary limit. In the case of our circularly symmetric 
water hole the stationary limit forms a ring. Bad for the boat: The stationary 
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limit is also the ring of no return. At best, the boat remains at its position, 
it never will escape. Any millimeter across the stationary limit will doom the 
boat, it will be inevitably sucked into the throat. Accordingly, the stationary 
limit coincides in this spherically symmetric case with the so-called event 
horizon. 

Event horizon 

In 1958, Finkelstein[ ] characterized the surface r = 2m as a “semi-permeable 
membrane” in spacetime, that is, a surface which can be crossed only in one 
direction. As soon as our boat has passed the event horizon, it can never 
come back. This property can be formulated in an invariant way: The light 
cones at each point of the surface have to nestle tangentially to the mem¬ 
brane. In 1964, Penrose [149] termed the null cone which divides observable 
from unobservable regions an event horizon. Mathematically speaking, the 
event horizon is characterized by having tangent vectors which are light-like 
or null at all points. Therefore, the event horizon is a null hypersurface. This 
is what is meant by a trapped surface[ 29], see Fig. 8 and Fig. 9, left image: a 
compact, spacelike, 2-dimensional submanifold with the property that out¬ 
going future-directed light rays converge in both directions everywhere on 
the submanifold. All these characterizations quite intuitively show up in the 
Penrose-Kruskal diagram to be discussed later. 

In view of the preceding paragraph, we define a black hole as a region 
of spacetime separated from infinity by an event horizon, see Carroll [29] and 
Brill[22]. 

Observational evidence in favor of black holes was reviewed by Narayan 
and McClintock[129]. 

Killing horizon 

The stationary limit surface is rendered more precise in the notion of a Killing 
horizon. A particle at rest (with respect to the infinity of an asymptotically 
flat, stationary spacetime) is to be required to follow the trajectories of the 
timelike Killing vector 9 . However, if we have a Killing vector K describing 

9 Using the definitions of the covariant derivative and of the Christoffel symbols, we 
can derive the following equation for an arbitrary vector K, 

K°d a g IJM = 2V Ul K v) - 2 g a ^d v) K a . (39) 

Assuming K a and g to be constant in time, demands = 0. Hence K has to be 

a Killing vector. In this coordinate system, we have K a K a = goo . Although K acts as 
time translation, it is not necessarily timelike! 
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a stationary spacetime, then at some points K may become lightlike, that 
is A' ,J /T /X = 0. If all these points build up a hypersurface E, then this null 
hypersurface is called a Killing horizon. Apparently, this notion is of a local 
character, in contrast to the definition of an event horizon, the definition of 
which refers to events in the future, it is of a nonlocal character, see Fig. 8. 

As we will see for the Schwarzschild black hole, see Fig. 9, outside the 
black hole the Killing vector is timelike, that is, K ll K fl < 0, on the Killing 
horizon it becomes null = 0 (by definition of the horizon), and inside 

it becomes spacelike A'''A), > 0. 

In the Schwarzschild case it will turn out that the event horizon and 
Killing horizon coincide, in the Kerr case they separate. 

Surface gravity 

From the definition of the Killing horizon it can be shown [29] that the quan¬ 
tity 

« 2 :=~(V a JC,)(VAA 1 ') | E (40) 

is constant on the Killing horizon and positive. The quantity k is called 
surface gravity. In simple cases, it has the interpretation of an acceleration 
or gravitational force per unit mass on the horizon. In the Schwarzschild 
spacetime it takes the value k = 1/4 m, which is the acceleration of a particle 
with unit mass as seen from infinity, compare with the Newtonian “field 
strength” (2) for r = 2m: 


GM m 1 
r 2 (2m) 2 4m 

In general, there is no such simple interpretation. 


(41) 


Infinite redshift 


Another property associated with the surface K 1 ' K u = 0 is the infinite red- 
shift. In view of the relation for the general relativistic time delay, 


r 0 {x B ) 


y/gtt(x B ) 
Vstt(x A ) 


r 0 (x A ). 


g tt —> 0 can be interpreted as follows. Consider tq(xb) the time measured by 
a clock B resting well away from the Killing horizon, whereas clock A with 
t 0 (x a ) is nearly at the Killing horizon. If gtt(x A ) —> 0 we get Tq(x b ) —> oo. 
From the point of view of clock B, clock A’s last signal, right before A hits 
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Figure 8: A null hypersurface is not necessarily an event horizon: Imagine 
a light cone that touches a hypersurface along the line of contact. Thus, the 
light cone is tangent as well as normal (in a spacetime sense) to the surface. 
Consequently, all such surfaces are null hypersurfaces. In the cases A and 

B, the light cone is entirely trapped inside the surface. Case A suggests that 
the surface does not close in a finite region, therefore the enclosed volume 
is not compact. Case B represents a (part of a) circle, which encloses all 
tangential light cones, and this forms an (black hole) event horizon. In case 

C, the light cone intersects the hypersurface. The white domain is outside 
the null surface but inside the light cone and, thus, reachable from within 
the enclosed domain. 
























the Killing horizon, will not reach B in a finite time, that is, never. To put it 
a little bit different: Signals sent with respect to A with constant frequency 
arrive increasingly delayed at B; for B the frequency approaches zero. This 
is called infinite redshift. 

Let us work out these ideas for the Schwarzschild solution and let us take 
“photons” in spacetime instead of boats on a lake. 


2.5 Using light rays as coordinate lines 

Schwarzschild coordinates exhibit a coordinate singularity at r = 2m. This 
obstructs the discussion of the event horizon considerably. As we have seen, 
light rays penetrate the horizon without difficulty. This suggests to use light 
rays as coordinate lines. Therefore we introduce in- and outgoing Eddington- 
Finkelstein coordinates. By combining both, we arrive at Kruskal-Szekeres 
coordinates, which provide a regular coordinate system for the whole Schwarz¬ 
schild spacetime. 


Eddington-Finkelstein coordinates 


In relativity, light rays, the quasi-classical trajectories of photons, are null 
geodesics. In special relativity, this is quite obvious, since in Minkowski space 
the geodesics are straight lines and “null” just means v = c. A more rigorous 
argument involves the solution of the Maxwell equations for the vacuum 
and the subsequent determination of the normals to the wave surface (rays) 
which turn out to be null geodesics. This remains valid in general relativity. 
Null geodesics can be easily obtained by integrating the equation 0 = ds. 
We find for the Schwarzschild metric, specializing to radial light rays with 
df = 0 = dd, 


t = ± 


+ 2m In 



+ const. 


(42) 


If we denote with r 0 the solution of the equation r + 2m In ^ — 1 = 0, we 
have for the t-coordinate of the light ray t(ro) =: v. Hence, if r = ro, we can 
use v to label light rays. In view of this, we introduce v and u 


r 


t + r + 2m In 

/ 

-1 

2m 

r 

t — r — 2m In 

-1 

2m 


(43) 

(44) 


Then ingoing null geodesics are described by v = const , outgoing ones by 
u = const, see Fig. 9. We define ingoing Eddington-Finkelstein coordinates 
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Figure 9: In- and outgoing Eddington-Finkelstein coordinates (where we 
introduce if with v = t' + r, u = t' — r). The arrows indicate the direction 
of the original Schwarzschild coordinate time (and thereby the direction of 
the Killing vector d t ). The left figure illustrates a black hole: All incoming 
photons traverse the event horizon and terminate in the singularity. The 
right figure illustrates a white hole: All outgoing photons emerge from the 
singularity, cross the horizon, and propagate out to infinity. 

by replacing the “Schwarzschild time” t by v. In these coordinates (v, r, 9, 0), 
the metric becomes 

2777 , \ 

1-) dv 2 + 2 dv dr + r 2 dfl 2 . (45) 

r J 

For radial null geodesics ds 2 = d9 = d<f> = 0, we find two solutions of (45), 
namely v = const and v = Am In |r/2m — 1| + 2r + const. The first one de¬ 
scribes infalling photons, i.e., t increases if r approaches 0. At r = 2m, there 
is no singular behavior any longer for incoming photons. Ingoing Eddington- 
Finkelstein coordinates are particular useful in order to describe the grav¬ 
itational collapse. Analogously, for outgoing null geodesics take (u, r, 9 , 0) 
as new coordinates. In these outgoing Eddington-Finkelstein coordinates the 
metric reads 

ds 2 = — ^1-^ du 2 — 2 du dr + r 2 dQ 2 . (46) 

Outgoing light rays are now described by u = const , ingoing light rays by 
u = —(4m In |r/2m — 1| + 2 r) +const. In these coordinates, the hypersurface 
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r = 2 m (the “horizon”) can be recognized as a null hypersurface (its normal 
is null or lightlike) and as a semi-permeable membrane. 


Kruskal-Szekeres coordinates 


Next we try to combine the advantages of in- and outgoing Eddington- 
Finkelstein coordinates in the hope to obtain a fully regular coordinate 
system of the Schwarzschild spacetime. Therefore we assume coordinates 
(u,v,9,4>)- Some (computer) algebra yields the corresponding representa¬ 
tion of the metric: 

ds 2 = — ( 1- 7 -- ] du dv + r 2 (u, v) dfl 2 . (47) 

V r{u,v)J 

Unfortunately, we still have a coordinate singularity at r = 2 m. We can get 
rid of it by reparametrizing the surfaces u = const and v = const via 

5 = exp (£)’ fi = - exp (-£) ■ (48) 

In these coordinates, the metric reads (r = r(u,v) is implicitly given by (48) 
and (44), (43), r s = 2m) 

di , 2 = --kLexp(-dtlt)) dvdu + r 2 (u,v)d(l 2 . (49) 

r(u, v) \ 2m ) 


Again, we go back from u and v to time- and space-like coordinates: 


t : = - (v + u) , 


r := 2 {v-u) 


In terms of the original Schwarzschild coordinates we have 

— 1 exp ( -— ) cosh 


r = 


t = 


2m 


2m 


- 1 


exp 

exp 


\4m7 

\4m/ 


4m ’ 


( 


4m 


(50) 


(51) 

(52) 


The Schwarzschild metric 

ds 2 = exp ( —(— dt 2 + dr 2 ) + r 2 dQ 2 , (53) 

r V 2m / v ' 

in these Kruskal-Szekeres coordinates behaves regularly at the 

gravitational radius r = 2m. If we substitute (53) into the Einstein equation 
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(via computer algebra), then we see that it is a solution for all r > 0. Eqs.(51) 
and (52) yield 


r — t' 




(54) 


Thus, the transformation is valid only for regions with |f| > t. However, we 
can find a set of transformations which cover the entire (t, r)-space. They are 
valid in different domains, indicated here by I, II, III, and IV, to be explained 
below: 


(I) 

(II) 
(HI) 

(IV) 


* = V^T^-p(^)sinh^ 

? = exp (£) cosh £ 

t: = ex p(i^) cosh £ 

f = V 71 -^ ex P ( 4 ^) sinh £ 

t = -V 2 H - 1 ex P (in) sinh £ 

' = “V 7 ^ -p (£) cosh V. 

1 = -v 71 -^ ex P (in) cosll £ 

? = -V^^exp^sinh^ 


The inverse transformation is given by 

(— -l) exp (——) 
\2 m ) F \2 m) 

f J Artanh t/r , 

Am 1 Artanh r/t, 


~2 72 

— r —t , 

for (I) and (III), 
for (II) and (IV) . 


(55) 

(56) 

(57) 

(58) 

(59) 

(60) 


The Kruskal-Szekeres coordinates (t, r, 9 , </>) cover the entire spacetime, see 
Fig. 10. By means of the transformation equations we recognize that we need 
two Schwarzschild coordinate systems in order to cover the same domain. 
Regions (I) and (III) both correspond each to an asymptotically flat universe 
with r > 2 771 . Regions (II) and (IV) represent two regions with r < 2m. Since 
t is a time coordinate, we see that the regions are time reversed with respect 
to each other. Within these regions, real physical singularities (corresponding 
to r = 0) occur along the curves t 2 — r 2 = 1. From the form of the metric 
we can infer that radial light-like geodesics (and therewith the light cones 
ds = 0) are lines with slope 1. This makes the discussion of the causal 
structure particularly simple. 
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Figure 10: Kruskal-Szekeres diagram of the Schwarzschild spacetime. 

2.6 Penrose-Kruskal diagram 

We represent the Schwarzschild spacetime in a manner analogous to the Pen¬ 
rose diagram of the Minkowski spacetime. To this end, we proceed along the 
same line as in the Minkowskian case. 

First, we switch again to null-coordinates v' = t + r and v! = t — f 
and perform a conformal transformation which maps infinity into the finite 
(again, by means of the tangent function). Finally we return to a time-like 
coordinate t and a space-like coordinate f. We perform these transformations 
all in one according to 

t + f — tan - , t — f = tan—— . (61,62) 
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Figure 11: Penrose-Kruskal diagram of the Schwarzschild spacetime. Region 
II corresponds to a black hole, region IV to a white hole. Regions I and III 
correspond to two universes. 


The Schwarzschild metric then reads 


„3 exp 
ds 2 = rs 


r(r,t) 
2 m 


(— dt 2 + dr 2 ) 


r(r,t ) 


cos 2 cos 2 Lp 


+ r 2 (i, f) dQ 2 , 


(63) 


where the function r(t, r) is implicitly given by 


( r \ ( r \ t + f t — r 

\2m ) f V2 m) 2 2 v ’ 

The corresponding Penrose-Kruskal diagram is displayed in Fig. 11 . The 
notations for the different infinities can be extracted from Table 1. In contrast 
to Minkowski space, light rays and particles may not escape to infinity but 
enter the black hole (II). Likewise, light rays and particles may not emerge 
from infinity but from the white hole (IV). 
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2.7 Adding electric charge and the cosmological con¬ 
stant: Reissner Nordstrom 


As mentioned in the historical remarks, soon after Schwarzschild’s solution, 
the first generalizations, including electric charge and the cosmological con¬ 
stant were published. We can be even quicker ... We already calculated the 
Ricci tensor for the Gullstrand-Painleve ansatz. If we use the well-known 
energy-momentum tensor for a point charge [81], the field equation may be 
written as 10 


V - + A y = 'tA 0 ^ i diag(-l,-1,1,1). 


(65) 


Taking the trace, we find R = 4A and arrive at 


P 2 _ p s _ 2 dr-NO _ q 2 

1X2 — ^3 — -o- — 2 A -7 • 


( 66 ) 


This equation can be integrated elementarily, 
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q 2 2a 
^ ~ 


(67) 


This function also solves the remaining two field equations. The integration 
constant a is again the mass m. Substituted into (37) and transformed to 
Schwarzschild coordinates (/ = 1 — 20) the solution reads 


ds 2 = —f(r ) dt 2 + + r 2 dQ 2 , 

fir) 


with 


ft \ i 2m , d 2 A 2 

/(r) := 1- + r 


( 68 ) 


(69) 


A detailed derivation using Schwarzschild coordinates and computer algebra 
can be found in Puntigam et al. [158] 

A discussion of the Reissner-Nordstrom(-de Sitter) solution can be found 
in Griffiths & Podolsky [69], for example. We only remark, that we recover 
the Schwarzschild solution for q = 0 and A = 0. The algebraic structure of 
the solution is identical to the Schwarzschild case. Thus, we find, in general, 
a singularity at r = 0. ffowever, a pure cosmological solution, m = 0, q = 0 


10 Einstein’s gravitational constant is denoted by n, Ao = y ^ is the admittance of the 
vacuum. With c = 1 and G = 1 we have «An = 2. 
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and A / 0, possesses no singularity and no horizon! On the other hand, an 
electrically charged black hole, A = 0, exhibits two horizons, 

f{r) = 0 <=>■ r± = m± \Jm 2 — q 2 . (70) 

In this respect, the charged black hole shows some similarities to a rotating 
(Kerr) black hole. We will pick up this discussion in Sec. 3.4. 


2.8 The interior Schwarzschild solution and the TOV 
equation 

In the last section we investigated the gravitational field outside a spherically 
symmetric mass-distribution. Now it is time to have a look inside matter, 
see Adler, Bazin, and Schiffer [ ]. Of course, in a first attempt, we have 
to make decisive simplifications on the internal structure of a star. We will 
consider cold catalyzed stellar material during the later phase of its evolution 
which can be reasonably approximated by a perfect fluid. The typical mass 
densities are in the range of ~ 10'g/cm 3 (white dwarfs) or « 10 14 g/cm 3 
(neutron stars, e.g., pulsars). In this context we assume vanishing angular 
momentum. 

We start again from a static and spherically symmetric metric 

ds 2 = -e^ (r) c 2 dt 2 + e B(r) dr 2 + r 2 dQ 2 (71) 


and the energy-momentum tensor 

Tiu, = (p + u„u u +pg^„, (72) 

where p = p(r) is the spherically symmetric mass density and p = p(r) the 
pressure (isotropic stress). This has to be supplemented by the equation of 
state which, for a simple fluid, has the form p = p(p). 

We compute the non-vanishing components of the field equation by means 
of computer algebra as (here k = 8^ rG/c 4 is Einstein’s gravitational constant 
and ()' = d/dr) 

— e B nr 2 c 2 p + e B + B'r — 1 = 0 , (73) 

—e B npr 2 — e B + A'r + 1 = 0 , (74) 

-Ae B npr+ 2A"r + {A’f r - A’B’r+ 2A 1 -2B' = 0. (75) 


The (0, 0)-component turns out to be equivalent to the (0, 0)-component. 
For convenience, we define a mass function m(r) according to 


-B 


=: 1 — 


2 m(r) 

r 


(76) 
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We can differentiate (76) with respect to r and find, after substituting 
(73), a differential equation for m(r) which can be integrated, provided p(r) 
is assumed to be known 


m(r ) = J | pc 2 £ 2 d£ . (77) 

Differentiating (74) and using all three components of the field equation, we 
obtain a differential equation for A: 


A' = 


2/7 

p + pc 2 


(78) 


We can derive an alternative representation of A 1 by substituting (76) into 
(74). Then, together with (78), we arrive at the Tolman-Oppenheimer-Volkoff 
(TOV) equation 

/ = (pc 2 +p)(m + hi pr 3 / 2) 

^ r(r — 2m) 

Terms that survive in the Newtonian limit are emphasized by boldface letters. 
The system of equations consisting of (77), (78), the TOV equation (79), 
and the equation of state p = p(p) forms a complete set of equations for the 
unknown functions A(r), p(r), p(r), and m(r), with 


dr 2 

ds 2 = — e A(r ^ c 2 dt 2 H-+ r 2 dVL 2 . (80) 

These differential equations have to be supplemented by initial conditions. 

In the center of the star, there is, of course, no enclosed mass. Hence we 
demand m(0) = 0. The density has to be finite at the origin, i.e. p( 0) = p c , 
where p c is the density of the central region. At the surface of the star, 
at r = Rq, we have to match matter with vacuum. In vacuum, there is 
no pressure which requires p(Rq ) = 0. Moreover, the mass function should 
then yield the total mass of the star, m(R & ) := GM & /c 2 . Finally, we have 
to match the components of the metric. Therefore, we have to demand 
exp[A(r 0 )] = 1 - 2 m(R e )/R e . 

Equations (73), (74), (75) and certain regularity conditions which gener¬ 
alize our boundary conditions, that is, 


• regularity of the geometry at the origin, 


• finiteness of central pressure and density, 

• positivity of central pressure and density, 
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positivity of pressure and density, 


• monotonic decrease of pressure and density, 

impose conditions on the functions p and p. Then, even without the explicit 
knowledge of the equation of state, the general form of the metric can be 
determined. For recent work, see Rahman and Visser [162] and the literature 
given there. 

We can obtain a simple solution, if we assume a constant mass density 

p — d( r ) = const. (81) 


One should mention here that p is not the physically observable fluid density, 
which results from an appropriate projection of the energy-momentum tensor 
into the reference frame of an observer. Thus, this model is not as unphysical 
as it may look at the first. However, there are serious but more subtle 
objections which we will not discuss further in this context. 

When p = const., we can explicitly write down the mass function (77), 


m(r = 


with R = 


m Q := 




2 R 2 V i'Y" ” 2 R 2 

This allows immediately to determine one metric function 

1 


e B = 


1 - 


R 2 


The TOV equation (79) factorizes according to 

| = -^(pc 2 + P )(i + K S‘ P ) ¥ L- 

It can be elementarily solved by separation of variables, 

Jfc-Rl-Vfc-r* 

pyr) = pc' —- =. 

Pa 1 - r 2 - 3 Jr 2 - 


(82) 


(83) 


(84) 


(85) 


Using (78) as A' = —2 [ln(p + pc 2 )}' and continuous matching to the exte¬ 
rior, eventually yields the interior & exterior Schwarzschild solution for a 
spherically symmetric body [173] 


ds 2 = 


(id 1 -#-Id 1 -!;) C 2 df- + j-y drC rW , r < R 0 
\ / r 2 


T 

r 2 


R 2 


1 

i 

2m0 


( 1 _ 2 ^ 0 ) c 2 dt 2 + _l_ dr 2 +r 2 dn 2^ 
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r>R e 

( 86 ) 


















The solution is only defined for R & < R. For the Sun 11 we have M 0 « 
2 x 10 30 kg, R q « 7x 10 8 m and accordingly p Q ta 1.4 x 10 3 kg/m 3 . This 
leads to R ~ 3 x 10 11 m, that is, the radius of the star i? 0 is much smaller 
than R: R & < R. Hence the square roots in (86) remain real. 

The condition R Q < R suggests that a sufficiently massive object cannot 
be stable since no static gravitational field seems possible. This conjecture 
can be further supported. Even before reaching R, the central pressure be¬ 
comes infinite, 



(87) 


If there is no static solution and the situation remains spherically symmetric, 
we are forced to the conclusion that such a mass distribution must radially 
collapse; either in an infinite time or to a single point in space. With reason¬ 
able simplifications, it was first shown by Oppenhcimer and Snyder[145] that 
the second alternative is true: A very massive object collapses to a black 
hole. As various singularity theorems show today, this behavior is indeed 
generic, see Chruscicl et al. [38] and Sec.3.10. 


u To ascertain the consistency of dimensions and units, we recollect the basic definitions: 


(m/s) 4 m 3 


8t tG 


N kgs 2 ’ 


The mass M carries the unit kg, the mass parameter has the dimension of a length: 


GM m 3 kp , s 2 



The definition of m(r) in Eq.(77) is consistent. For p = const , we have 



GM(r ) 


Here p denotes the physical mass density, \p] = kg/m 3 . Thus 


4 

M(r) := -7 t r 6 p 

O 

is the physical mass with the unit kg. 
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3 The Kerr metric (1963) 

After some historical reminiscences (Sec.3.1), we point out how one can ar¬ 
rive at the Kerr metric (Sec. 3.2). For that purpose, we derive, in cylindrical 
coordinates, the four corresponding partial differential equations and explain 
how this procedure leads to the Kerr metric. In Sec. 3.3, we display the Kerr 
metric in three classical coordinate systems. Thereafter we develop the con¬ 
cept of the Kerr black hole (Sec. 3.f). In Secs.3.5 to 3.1, we depict and dis¬ 
cuss the geometrical/kinematical properties of the Kerr metric. Subsequently, 
in Sec. 3.8, we turn to the multipole moments of the mass and the angular 
momentum of the Kerr metric, stressing analogies to electromagnetism. In 
Sec.3.9, we present the Kerr-Newman solution with electric charge. Even¬ 
tually, in Sec.3.10, we wonder in which sense the Kerr black hole is distin¬ 
guished from the other stationary axially symmetric vacuum spacetimes, and, 
in Sec.3.11, we mention the rotating disk metric of Neugebauer-Meinel as a 
relevant interior solution with matter. 


....When I turned to Alfred Schild, who was still sitting in 
the armchair smoking away, and said “Its rotating!” he was 
even more excited than I was. I do not remember how we 
celebrated, but celebrate we did! 

Roy P. Kerr (2009) 


3.1 Historical remarks 

The search for axially symmetric solutions of the Einstein equation started 
in 1917 with static and was extended in 1924 1° stationary metrics. It cul¬ 
minated in 1963 with the discovery of the Kerr metric. 

The Schwarzschild solution, as we have seen, describes the gravitational 
field of a spherically symmetric body. Obviously, most planets, moons, and 
stars rotate so that spherical symmetry is lost and one spatial direction is 
distinguished by the 3-dimensional angular momentum vector J of the body. 
Hence the next problem to attack was to search for the gravitational field of 
a massive rotating body. 

When one considers a static and axially symmetric situation—this is the 
case if the body does not carry angular momentum—then one can choose the 
rotation axis as the 2 -axis of a cylindrical polar coordinate system: x l = z, 
x 2 = p and x 3 = 0. Then static axial symmetry means that the components 
of the metric g^ u = g^ u (z,p) do not depend on the time t and the azimuthal 
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angle <ft (we have here one timelike and one spacelike Killing vector 12 ). 

Already in 1917, Weyl[188] started to investigate static axially symmetric 
vacuum solutions of Einstein’s field equation. He took cylindrical coordinates 
and proposed the following “canonical” form of the static axisymmetric vac¬ 
uum line element: 13 


ds 2 = fdt 2 — |/r(d; 2 2 + dp 2 ) + ~E~ _ | > (89) 

here / = f(z,p ) and h = h(z,p ) and (x° = t, x 1 = z, x 2 = p, x 3 = (ft). 
Weyl was led, in analogy to Newton’s theory, to a Poisson equation and 
found thereby a family of static cylindrically symmetric solutions that could 
be understood as the exterior field of a line distribution of mass along the 
rotation axis. Similar investigations were undertaken by Levi-Civita[108l 
(1917/19). 

In the year 1918, Lense and Thirring[107] investigated a rotating body. 
They specified the energy-momentum tensor of a slowly rotating ball of mat¬ 
ter of homogeneous density and integrated the Einstein equation in lowest 
approximation. They found, for a ball rotating around the z-axis of a spa¬ 
tial Cartesian coordinate system, the linearized Schwarzschild solution in 
isotropic coordinates, see Table 2, together with two new “gravitomagnetic” 
correction terms in off-diagonal components of the metric (k is Einstein’s 
gravitational constant): 

I (dx 2 + dy 2 + dz 2 ) - -^(xdy — ydx)dt 

> — r ;___- 

linearized Schwarzschild gravitomagnetic term 

(90) 

12 Remark on Killing vectors: Consider a point P of spacetime with coordinates x a . 
We specify a direction at P. If we have a flat Minkowski space, the components g^ v 
of the metric, given in Cartesian coordinates, would not change under a motion in the 
^-direction. However, in a curved spacetime, the g^ will change in general. If fulfills 
the Killing equations (see Stephani [179]) 

= 0 , ( 88 ) 

with V as covariant derivative operator, then is called a Killing vector , and this vector 
specifies a direction under which the metric does not change. The Schwarzschild metric is 
static, that is, it has one timelike Killing vector along the time coordinate. Furthermore, 
it is spherically symmetric and thus has three additional spacelike Killing vectors. In 
the Weyl case, because of the axial symmetry around the z-axis, two of those spacelike 
Killing vectors get lost. Left over in the Weyl case are the two Killing vectors, one timelike 
= Q t anc [ one spacelike = d^. 

13 Weyl used p —► r ,</>—» •&. 
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this is valid for kM <C r and kJ z <C r 2 . This gravitomagnetic effect (“the 
Lense-Thirring effect”) is typical for GR: in Newton’s theory a rotating rigid 
ball has the same gravitational field as a non-rotating one. Gravitomagnetism 
is alien to Newton’s gravitational theory. 

In the meantime, the Lense-Thirring effect has been experimentally con¬ 
firmed by the Gravity Probe B experiment, see Everitt et al.[ ]. They took 
a gyroscope in a satellite falling freely around the (rotating) Earth. The spin 
axis of the gyroscope pointed to a fixed guide star. Because of the gravito¬ 
magnetic term in (90), the gyroscope executed a (very small) Lense-Thirring 
precession. 11 This can be understood as an interaction of the spin of the 
gyroscope with the spin of the Earth (spin-spin interaction). Since the gy¬ 
roscope moves along a 4d geodesic of a spacetime curved by the mass of the 
Earth, an additional geodetic precession occurs that has to be experimen¬ 
tally separated from the Lense-Thirring term. The geodetic precession had 
already been derived earlier by de Sitter[ I ] in 1916. 15 

In spherical polar coordinates we have ydx — xdy = r 2 sin 2 # d(p. Thus, the 
gravitomagnetic cross-term in (90) may be rewritten as (AkJ z s\\\ 2 6/r) dt d(j). 
A comparison with (89) shows that the canonical Weyl form of the static 
metric is too narrow for describing rotating bodies. 

From 1919 on, there appeared further articles on axisymmetric solutions. 
Levi-Civita 16 (1919) reacted to Weyl’s article, and Bach[8] (1922) pushed the 
Lense-Thirring line element to the second order in the approximation. 

Then, in 1924, Lanczos[105] extending the Weyl ansatz, started to inves¬ 
tigate stationary 17 solutions. He found an exact solution for uniformly ro¬ 
tating dust. However, his work was apparently partially overlooked. Later, 
Akeley[2, 3] (1931), Andress[6] (1930) and, in a more definite form, Lewis[109] 
(1932) generalized the static Weyl metric to a stationary one by taking into 
account the gravitomagnetic term of Lense-Thirring. Lewis (1932) wrote, in 
cylindrical polar coordinates (x\ p , x 2 z), 

ds 2 = fdt 2 — (e^dxi 2 + e u dx 2 2 + ld(j) 2 ) — 2 mdt d<f>. (91) 

He found some exact solutions, typically for rotating cylinders, but not for 

14 For related experiments, see Ciufolini et al.[35, 34] and Iorio et al.[89, 92] A recent 
comprehensive review was given by Will[189]. A textbook presentation may be found in 
Ohanian & Ruffini[143]. 

15 De Sitter had applied it to the Earth-Moon system conceived as a gyroscope precessing 
around the Sun (the rotation of which can be neglected). This effect can nowadays be 
measured by Lunar Laser Ranging, see Will[L 89] . 

16 See Ref. [108], note 8 with the subtitle “Soluzioni binarie di Weyl”. 

17 Stationary spacetimes are those that admit a time-like Killing vector. Static space- 
times are stationary spacetimes for which this Killing vector is hypersurface orthogonal; 
physically this implies time reversal invariance and thus the absence of rotation. 
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rotating balls. It became definitely clear that, in the axially symmetric 
case, we may have many different exact vacuum solutions, in contrast to 
the case of spherical symmetry with, according to the Birkhoff theorem, the 
Schwarzschild solution as being unique. 

Not much later, van Stockum[184] (1937) determined the gravitational 
field of an infinite rotating cylinder of dust particles, thereby recovering the 
Lanczos solution, inter alia. He fitted one of the interior matter solutions of 
Lewis to an exterior vacuum solution. Continuing on this line of research, 
Papapetrou[148] (1953) started from the Andress-Lewis line element, putting 
it in a slightly different form, suitable for all stationary axisymmetric vacuum 
solutions: 


ds 2 = —e^(dp 2 + dz 2 ) — ld(j) 2 — 2m dcj) dt + fdt 2 . (92) 

The functions //, l, m, and / depend only on p and z. Papapetrou integrated 
the field equations and found exact stationary rotating vacuum solutions. 
However, his solution carried either mass and no angular momentum or an¬ 
gular momentum and no mass. Thus[148], “this solution is very special and 
physically of little interest.” 

A year later, a new result was published, which gave the problem of 
finding solutions for a rotating ball a new direction. Petrov[152] (1954), 
from Kazan, classified algebraically the Einstein vacuum field, that is, the 
Weyl curvature tensor, according to its eigenvalues and eigenvectors. This 
information reached the West, in the time of the Cold War, with some delay. 
A bit later, Pirani[154] (1957) developed a related formalism. It was the 
Petrov classification and the picking of a suitable class for the gravitational 
field of an isolated body (Petrov class D, with two double principal null 
directions) that finally led to the discovery of the Kerr solution during 1963, 
ten years after the unphysical solutions of Papapetrou. 

Accordingly, it turned out to be a formidable task to find an exact solution 
for a rotating ball and it was only found nearly half a century after the 
publication of Einstein’s field equation, namely in 1963 by Roy Kerr [98], a 
New Zealander, who worked at the time in Texas within the research group of 
Alfred Schild. It is a 2-parameter solution of Einstein’s vacuum field equation 
with mass M and rotation (or angular momentum) parameter a := J/M. 

The story of the discovery of the Kerr solution was told by Kerr himself 
at a conference on the occasion of his 70th birthday [99]. A decisive starting 
point of Kerr’s investigations was, as mentioned, the Petrov classification. 
Melia, in his popular book[L24] “Cracking the Einstein Code”, which does 
not contain any mathematical formula—apart from those appearing in two 
copies of Kerr’s notes and on a blackboard in another figure—has told this 
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fascinating battle for solving Einstein’s equation, see also the Kerr story in 
Ferreira [58]. 

Dautcourt [10] discussed the work of people who were involved in this 
search for axially symmetric solutions but who were not so fortunate as 
Kerr. In particular, Dautcourt himself got this problem handed over from 
Papapetrou in 1959 as a subject for investigation. He used the results of Pa- 
papetrou (1953). Dautcourt’s scholarly article is an interesting complement 
to Melia’s book. In particular, it becomes clear that the (Lanczos-Akeley- 
Andress-Lewis-)Papapetrou line element (92) was the correct ansatz for the 
stationary axially symmetric metric and the Kerr metric is a special case 
therefrom. The Papapetrou approach with the line element (92) was later, 
after Kerr’s discovery, brought to fruition by Ernst [53] and by Kramer and 
Neugebauer[100]. 

3.2 Approaching the Kerr metric 

We derive a 2nd order partial differential equation, the Ernst equation, that 
governs the stationary axially symmetric metrics in Einstein’s theory. Subse¬ 
quently, we sketch how the Kerr solution emerges as a simple case therefrom. 

Papapetrou line element and vacuum field equation 

In more modern literature, the Papapetrou line element (92), which describes 
some rotation around the axis with p = 0, is usually parametrized as fol¬ 
lows 18 , 

ds 2 = f(dt — udf) 2 — / _1 \e 2l (dp 2 + dz 2 ) + p 2 d(j) 2 ] , (93) 

te (—00,00), p e [ 0 ,oo), z e (— 00,00), <f>&[ o, 27 t); 

we assume / > 0. We compute the vacuum field equation of this metric. 
Nowadays we can do this straightforwardly with the assistance of a computer 
algebra system. During the 1960s, when this work was mainly done, there 
were no computer algebra systems around. Hearn[79] released the computer 
algebra system REDUCE in 1968. Back then, one had to be in command of 
huge computer resources in order to bring the underlying computer language 
LISP to work. Today, Reduce can run on every laptop; for other computer 
algebra systems, see Grabmeier et al. [68] and Wolfram[191], 

Because of its efficiency, we will use Schriifer’s Reduce-package EXCALC, 
which was built for manipulating expressions within the calculus of exterior 

18 See Ernst [53], Buchdahl[23], de Felice & Clarke[57], Quevedo [160], 0’Neill[ 1], 
Stephani et al. [180] , Eq.(19.21), Griffiths & Podolsky[69], and Sternberg[181]. 


44 



forms. For that purpose, we reformulate the metric (93) in terms of an or¬ 
thonormal coframe of four 1 -forms d a = e* a dx\ with the unknown functions 
/ = f{p,z), u = u(p,z), and 7 = 7 (p,z), namely 


d® 

= / 5 (dt — ood 4>) 

= e^dx 1 — ft (dx° — udx 3 ), 

(94) 

d 1 

= f~^e 1 dp 

= e}dx l = f~^e 1 dx 1 , 

(95) 

d 2 

= f~^e y dz 

= ei 2 dx l = f~^e 1 dx 2 , 

(96) 

d 3 

= f~bd<l> 

= e^dx 1 = / _ 2 pdx 3 . 

(97) 


Because of the orthonormality of the coframe d a , we have 

ds 2 = g = + d° (g) d° - d 1 ® tf 1 - d 2 ® d 2 - d 3 ® d 3 . (98) 

Eqs.(94) to (98) are equivalent to (93). The corresponding computer code, 
as input for Reduce-Excalc, reads as follows: 
pform f=0, omega=0, gamma=0 $ 

fdomain f=f(rho,z), omega=omega(rho,z), gamma=gamma(rho,z); 


coframe o(0) = sqrt(f) * (d t - omega * d phi), 

o(l) = sqrt(f)**(-1) * exp(gamma) * d rho, 
o(2) = sqrt(f)**(-l) * exp(gamma) * d z, 
o(3) = sqrt(f)**(-1) * rho * d phi 
with signature (1,-1,-1,-1); 

Isn’t that simple enough? From this data, the Einstein equation is cal¬ 
culated, with the Einstein tensor G^ v . The complete, fairly trivial pro¬ 
gram is documented in the Appendix. Note in particular that we used a 
ETf^X interface allowing us to output the expressions directly in BTgX. This 
computer output—without changing anything of the formulas—after some 
post-editing for display purposes, reads as follows: 


G° 0 := (4 • dp,pf ■ f ■ p 2 - 5 • d p f 2 ■ p 2 + 4 • d p f ■ f ■ p + 4 • d z , z f ■ f ■ p 2 

-5 • d z f 2 ■ p 2 - 4 ■ d P)P 7 • f 2 ■ p 2 - 4 • <9^7 ■ f 2 ■ p 2 + 3 • d p u 2 • / 4 


+3 • d z uj 2 • 

f 4 ) / (4 • e 2 " 

7 -/-p 2 ) , 





(99) 

/ •(2 • d p f 

• dpLd ■ p + 2 

• d z f ■ 3 z uj ■ 

P + dp, p 

uj- f ■ p- dpu 

•/ 


+d Z)Z u ■ f • 

p)/( 2-e 2 ' 7 

V), 





(100) 

0 dpf 2 • P 2 - 

■ dzf 2 ■ p 2 ~ 

4 • dp'y ■ f 2 

■ P 

- dpJ 2 ■ f 4 + <9 2 w 2 

•/ 4 ) 


/ (4 ■ e 2 ' 7 ■ j 

■•P 2 ), 






(101) 

( d p f ■ d z f ■ 

P 2 - 2 • 5^7 

■f 2 -p~d\ 

• 

d 2 w • 

/ 4 )/(2 • e 2 ' 7 • 

/•p 2 ; 

)(102) 

(- d P f 2 ■ p 2 

- d z f 2 -p 2 

- 4 • dp,pi ■ 

f 2 

•P 2 - 

- 4 • <9 2j2 7 • / 2 • 

P 2 



-d p u 2 • / 4 - d z u 2 • / 4 )/(4 • e 2 ' 7 • / • p 2 ). (103) 
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This calculation of the Einstein tensor by machine did not require more than 
about 15 minutes, including the programming and the typing in; for sample 
programs, see Socorro et al. [176] and Stauffer et al. [178]. 

Inspecting these equations, it becomes immediately clear that the numer¬ 
ator of (100) does not depend on 7 . In order to get a better overview, we 
abbreviate the partial derivatives of Reduce d p f by subscripts, f p , and drop 
the superfluous multiplication dots of Reduce. We find 

G 3 0 = 0 — >■ 0 = f(uj pp + uJ Z z — ~ u p) + 2 (/ p o;p + f z io z ). (104) 

Moreover, by subtracting (103) from (99) we find another equation free of 7 : 

G°o — G 3 3 = 0 — > 0 = f(f pp + -fp + fzz) ~ fp ~ f z + ^-(k’p+k’z) • (105) 

Left over are the equations (101) and (102), which can be resolved with 
respect to the first derivatives of 7 : 

G 1 i = 0 —>• 7 P = ~ fz ) + 

G 1 2 = 0 —> 7z = Ypfpfz ~ 2 ^ ) UJpUJz ' (107) 

Collected, we have these four equations determining the stationary ax- 
isymmetric vacuum metric: 


0 

— fifpp + - fp + fzz ) — fp — fz + “2 ( W P + w z) > 

(108) 

0 

fi^^pp Wzz ] p) “1“ 2 (^fpUJp + f z UJ z ) ? 

(109) 

7 p 

= ~ fz) + J~ p ^l - W p) , 

(110) 

7z 

p f f f 2 

- 2 pfpJz 2p UpUz- 

(111) 


Let us underline how effortless—under computer assistance—we arrived at 
these four partial differential equations (PDEs) for determining stationary 
axially symmetric solutions of Einstein’s field equation. 

Ernst equation (1968) 

It is one step ahead, before we arrive at a still more convincing form of these 
PDEs. After some attempts, one recognizes that (109) can be written as 

(W/(Wa 0 ' (112) 
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(113) 


With the ansatz (Q = h2(p, £)), 


f2 y 2 

^ p i ^ p 1 

p p 


f (fpp + -fp + fzz ) - fp - fz + - 0 • 


(114) 


Eq.(112) is identically fulfilled. We substitute (113) into (108): 

1 

P‘ 

Since (109) is already exploited, we can find fl by differentiating the f2’s 
in (113) with respect to 2 and p, respectively, and by adding the emergent 
expressions ( uj pz = co zp ): 


f (Qpp + ~^p + ^zz) ~ 2 fpQp ~ Zfz^z — 0. 


(115) 


Eqs.(108) and (115) can be put straightforwardly into a vector analytical 
form, if we recall that our functions do not depend on the angle 0: 19 


/A/-(V/)-V/ + (VU)-VU = 0, (116) 

/Aft-2(V/)-Vft = 0. (117) 


The last equation can also be written as V • (/~ 2 V Tl) = 0. Eqs.(116) and 
(117) liberate ourselves from the cylindrical coordinates, that is, this expres¬ 
sion is now put in form independent of the specific 3d coordinates. With the 
potential [i 2 = — 1 ) 

S:=f + iQ, (118) 


which was found by Ernst [53] and Kramer & Neugebauer[100], we find the 
Ernst equation[ >3] 

(Ref) AS = 'VS • 'VS , (119) 

or, in components, 


(Ref) 


~ d 2 £ 

dz 2 


ld_( d£Y 
P dp V dp)_ 



( 120 ) 


The “Re” denotes the real part of a complex quantity. Under stationary 
axial symmetry—the corresponding metric is displayed in (93)— the Ernst 
equation (119), together with Eqs.(118, 113, 110, 111), are equivalent to the 
vacuum Einstein field equation. 

19 In cylindrical coordinates, we have for a vector V and a scalar s the following formulas, 
see Jackson [94]: 


V-V = ±d p 0 oVP+dzVi + faVz, V 2 s 
Vs = ei d p s + e 2 d z s + ea^d^s, Vs-Vs 


As = l d p(P d P s ) + d 2 s+ jzdls, 
(< d p s ) 2 + (chs) 2 + ^(d^s) 2 . 
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From Ernst back to Kerr 

This reduces the problem of axial symmetry to the solution of the second 
order PDE (119). This method, which came along only five years after Kerr’s 
publication, led to many new exact solutions, amongst them the Kerr solution 
(1963) as one of the simplest cases. We are only going to sketch how one 
arrives at the Kerr solution eventually. We follow here closely Buchdahl [ ]. 
One introduces a new complex potential £ by 

f=; frr (121) 

Then the Ernst equation becomes 

(££-l)A£ = 2£V£-V£, (122) 


where the overline denotes complex conjugation. If one has a solution of this 
equation, we can determine the functions /, u and 7 by 


/ 

Up 


7 p 


Im[(£ + l) 2 £z] 

- l ) 2 


p --= 


Uz 


(£ - l ) 5 


UJ Z 


Im[(£ + 1) 2 £ P ] 

- l ) 2 


lz 


M.Uz) 

(« - 1) 2 ' 


(123) 

(124) 

(125) 


For rotating bodies, spherical prolate coordinates x, y, with a constant k, 
are much more adapted: 

p = k(x 2 — 1)^(1 — y 2 )^ , z — kxy. (126) 


It turns out that one simple potential solving the Ernst equation, with the 
constants p and q, is 


£ = px — iqy with p 2 + q 2 = 1; (127) 

it leads to the Kerr metric. For this purpose, one has to introduce the 
redefined constants m := k/p (mass) and a := kq/p (angular momentum 
per mass) and to execute subsequently the transformations px = (p/m) — 1 
and qy = (a/m) cos 9 to the new coordinates p and 6 . Then one arrives at 
the Kerr metric in Boyer-Lindquist coordinates, which is displayed in the 
table on the next page. For more detail, compare, for instance, the books of 
Buchdahl[23], Islam[93], Heusler[85], Meincl et al. [123], or Griffiths et al. [69] . 
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By similar techniques, a Kerr solution with a topological defect was found 
by Bergamini et al.[ ]. 

Incidentally, in the context of the Ernst equation, Geroch made the fol¬ 
lowing interesting conjecture: A subset of all stationary axially symmetric 
vacuum space-times, including all of its asymptotically flat members, that 
is, in particular the Kerr solution, can be obtained from Minkowski space by 
transformations generated by an infinite-dimensional Lie group. This con¬ 
jecture was “proved” by Hauser and Ernst [76], see also Ref. [ ] However, the 

proof contained a mistake that was subsequently corrected in Ref. [ ] 

Starting from 4d ellipsoidal coordinates, Dadhich[39] gave a heuristic 
derivation of the Kerr metric by requiring, amongst other things, that light 
propagation should be influenced by gravity. 

3.3 Three classical representations of the Kerr metric 

We collected these three classical versions of the Kerr metric in Table 4, see 
also Visser [ ] .Three more coordinate systems should at least be mentioned: 

• Pretorius & Israel[\ } double null coordinates: 
very convenient to tackle the initial value problem 

• Doran[ ] coordinates: 

Gullstrand-Painleve like; useful in analog gravity 

• Debever/Plebanski/Demianski[ 155] coordinates: 

components of the metric are rational polynomials; convenient for (com¬ 
puter assisted) calculations 
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Table 4. Kerr metric: the three classical representations 


Kerr-Schild 


(' t,x,y,z ) 


Cartesian background 


ds 2 = — dt 2 + dx 2 + dy 2 + cfe 2 


2 mr 3 / , r(xdx + ydy) a(ydx-xdy) z , 

-^ dt, + - n ~ + - cfe 


?’ 4 + a 2 z 2 


a z + r 2 


a 2 + r 2 


2,2,2 2,2(1 - I t 

x +y + z = r + a 1 -, , r = r(x, y, z 


(r cos 4 > + a sin </>) sin ( 
(r sin 4 > — a cos <f >) sin ( 
r cos 6 


Boyer-Lindquist 


(t,r,9,(j>) 


Schwarzschild like 


, n / 2mr\ ,9 4mrasin 2 # , ,, 

ds 2 = - 1 - dt 2 -=- dt d<f) 

V P ) P 

P 2 .2 , „2, fl 2 , (2 , „ 2 , 2 mra 2 sin 2 0 \ 2 


+ — dr 2 + p z d9 z + I r 2 + a z + 


sin 2 0 d(j) z 


p z := r 2 + a 2 cos 2 6 A := r z — 2mr + a 2 = (r — r-_)_)(r — ?’_) 


Kerr original 


r 2 + a 2 

an = dt H-—— ar 

dip = d<f> + — dr 


( v,r,e,tp ) 


Eddington-Finkelstein like 


ds 2 = - (l - (dv - a sin 2 9 dp) 2 


~\~ 2 (tfo — a sin 2 6 dcp'j (dv — a sin 2 9 dp) + p 2 ( d6 2 + sin 2 9 dp 2 ) 


rv , := m ± v m 2 — a 2 cos 2 0 r+ := m ± v ?n 2 — a 2 




As input for checking the Kerr solution, we use the orthonormal coframe[181] 

:= 

'd l 

d 2 := 


^ {dt — a sin 2 0 dp) , 

(128) 

P dr 

VS ’ 

(129) 

pdO , 

(130) 

-[(r 2 + a 2 ) dp — adt\ . 

(131) 


p 

We introduced the sign function, which is convenient for discussing the dif¬ 
ferent regions in the Penrose-Carter diagram: 


e = 


+1 for r > r + or r < r_ ,, 
-1 for r_ < r < r_i_. 


(132) 


The metric can then be written in terms of the coframe as 

ds 2 = g = e(-d° ® d° + d 1 <g> t? 1 ) + d 2 ® d 2 + d 3 ® d 3 . (133) 


From Table 4 it is not complicated to read off the Schwarzschild and the 
Lense-Thirring metric as special cases. In comparison to the Schwarzschild 
metric, the Kerr solution includes a new parameter a which will be related to 
the angular momentum. However, it should be noted that, by setting a = 0, 
the Kerr metric reduces to the Schwarzschild metric, as it should be (p 2 —> r 2 
and A —>• r 2 — 2 mr): 

ds 2 = — (1 — dt 2 H— tt— - dr 2 + r 2 d0 2 + r 2 sm 2 0d(j) 2 .(134) 

V r z ) r z — 2mr 


By canceling r 2 in the dr 2 - term, we immediately recognize the Schwarzschild 
metric. 

Considering the parameter a we should note the following fact. For small 
values of the parameter a, where we may neglect terms of the order of a 2 , we 
arrive at (p 2 —y r 2 and A —y r 2 — 2 mr) and 


ds 2 = 



dt 2 -|- dr 2 + r 2 (d0 2 + sin 2 6> dcp 2 ) 

r 


Ama sin 2 0 
r 


dt dtp. 


(135) 


Since, in spherical coordinates we have y dx—x dy = r 2 sin 2 8 dp, the crossterm 
may be rewritten as 4ma ^ in ~ e dtdp = (xdz — y dx). Thus, in the limiting 

case a 2 1, the Kerr metric yields the Lense-Thirring metric, provided we 
identify ma = J z . 
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3.4 The concept of a Kerr black hole 

We come back to our Fig. 7 with “Schwarzschild” versus “Kerr”. The Kerr 
spacetime may be visualized by a vortex, where the water of the lake spirals 
towards the center. Much of the above said for the Schwarzschild case is still 
valid. However, one important difference occurs. The stationary limit and 
the event horizon separate, which will be illustrated by corresponding graphical 
representations. 

In case of a vortex, the flow velocity of the in-spiraling water has two 
components. The radial component which drags the boat towards the center 
whereas the additional angular component forces the boat to circle around 
the center. Again, the stationary limit is defined by the distance at which 
the boat ultimately can withstand the radial and circular drag of the water 
flow. Beyond the stationary limit the situation is not as hopeless as in the 
Schwarzschild case. Using all its power, the boat may brave the inward flow. 
But then it has not enough power to overcome the angular drag and is forced 
to orbit the center. By means of a clever spiral course the boat may even 
escape beyond the stationary limit. The stationary limit is not necessary an 
event horizon. At some distance, nearer to the center than the stationary 
limit, also the pure radial flow of water will exceed the power of the boat. 
There, inside the stationary limit, is the event horizon. 

In order to investigate the structure of the Kerr spacetime, we first look at 
“strange behavior” of the metric components in Boyer-Lindquist coordinates. 
The following cases can be distinguished: 

A = 0 g rr becomes singular, 
p 2 = 2 mr g tt vanishes, 

p 2 = 0 g rr and goo vanish, the other components are singular. 

As we have extensively discussed in the previous section, singularities 
of components of the metric may signify physical effects but, on the other 
hand, may only be due to “defective” coordinates. Thus, we will proceed 
along similar lines to investigate the nature of these singularities. 

We will not address the geodesics of the Kerr metric in detail. For an 
elementary discussion the reader is referred to Frolov and Novikov[64] and 
to the more advanced discussion in Hackmann et al. [72]. 
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Depicting Kerr geometry 


We draw a picture of the spatial appearances and relations of the various 
horizons and the singularity of the Kerr metric. From outside to inside these 
are, explicitly, 


outer 

ergosurface 

r E+ 

:= m 

+ 

\Zm 2 — a 2 cos 2 9 

event 

horizon 

r+ 

:= m 

+ 

^ joined at polar axis 

\J m 2 — a 2 

Cauchy 

horizon 

r _ 

:= m 

— 

^ merge for a —> m 

\J m 2 — a 2 

inner 

ergosurface 

r E - 

\= m 

— 

^ joined at polar axis 

\/m 2 — a 2 cos 2 9 


singularity 

r 

= 0 


^ lies on the rim for 9 = 7t/2 


For a = 0, inner ergosurface and Cauchy horizon vanish, whereas outer 
ergosurface and event horizon merge to the Schwarzschild horizon. To vi¬ 
sualize the various surfaces we use Kerr-Schild quasi-Cartesian coordinates. 
The radial coordinate r of the Boyer-Lindquist coordinates is related to the 
coordinates x, y, z of the Kerr-Schild coordinates via, see Table 4, 

x + y -\ --— z = r + a , z = r cos 9 . (137) 

Substituting r = 0, r = r±, r = r E± , and a little bit of algebra yields: 

• Singularity r = 0 

Since r — 0 leads to z — 0, we get the equation of a circle of radius a 
in the equatorial plane, 

x 2 + y 2 = a 2 . (138) 

For a = 0, the ring collapses to the Schwarzschild singularity. 

A closer inspection shows that the structure of the singularity is more 
complex[66, 110]. 


Horizons r = r± 

In this case we arrive at the equation for an oblate (for a < m) ellipsoid, 


x 2 y 2 z 2 

—2 H-9 H-9 = 1 , 

ai «3 


(139) 


where a\ = a 2 = r± + a 2 > a 2 = \ . 
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a<m 



a=m 


event horizon 


ring singularity 



outer ergosurface 


ergoregion 


inner ergosurface 


ring singularity 



a>m 


(outer) ergosurface 


Fig. 12: Ergosurfaces, horizons, and singularity for slow, extremal (‘critical’), 
and fast Kerr black holes. 
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• Ergosurfaces r = r E ± (0) 

Things are a little bit more involved in this case because r is not con¬ 
stant. We can also derive a “ellipsoid-like” equation (for a < m), 

x 2 U~ z 2 

-1— -1-= 1 (140) 

am + a|(0) + 4(9) ’ 1 J 

now with 

«?(») = «?(») = r 2 E± W + a\ a&e) =-±-. (141) 

r E±\ U ) 

The 0-dependence will deform the ellipsoid. On the equatorial plane 
with 9 = 7T we have r E = 0. Hence, a x = a 2 = a and a 3 diverges. This 
results in a non-regular rim on which the ring singularity is located. 

For a > m, the rg ± is partly not defined, since the term under the 
square-root changes sign, and becomes negative if 

171 

cos (9) = — . (142) 

a 

This defines two rings with 9i = arccos (m/a), 0 2 = n — 9 1 and r = m. 
As a consequence, the outer ergosurface only extends to these rings from 
the outside, and the inner ergosurface up to the rings from the inside. 
The outcome is a kind of torus. The center-facing side is constituted by 
a part of the inner ergosurface (along with the ring singularity), whereas 
the outside facing parts are given by a part of the outer ergosurface. 

An extensive discussion of the embedding of the ergosurfaces into Eu¬ 
clidean space, together with corresponding Mathematica-programs, can 
be found in Ref. [114]. 

The surfaces are visualized in Fig. 12. We did not use a faithful embedding 
but rescaled axes in order to achieve better visibility. 

The presence of the term y/m 2 — a 2 requires the distinction of three differ¬ 
ent cases dependent on the values of the mass parameter m and the angular 
momentum parameter a: 

m > a slow rotation m = a critical rotation m < a fast rotation 

The slow rotating case shows the richest structure. Both ergosurfaces and 
both horizons are present and distinct from each other. As a approaches 
77i, the event and the Cauchy horizon draw nearer and nearer. At critical 
rotation, a = m, both horizons merge into one single event horizon with 
r = 77i. Eventually, for fast rotation a > m, the event horizon disappears 
and reveals the naked ring singularity which now is located at the inner edge 
of the now toroidal shaped (outer) ergosurface. 
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3.5 The ergoregion 

We explore the region between the outer ergosurface and the event horizon. 
There it is not possible to stand still, anything has to rotate, even the event 
horizon. The compulsory rotation in the ergoregion allows one to extract 
energy from the black hole. This so-called Penrose process leads to black hole 
thermodynamics. 


Constrained rotation 

The outer ergosurface, r = Re + , is defined by the equation goo(r E+ ) = 0. 
Thus, it is a surface of infinite redshift and a Killing horizon. For a third 
characterization of the ergo surfaces we have to deal not only with radial but 
also with rotational motion. Consider the Kerr metric in Boyer-Lindquist 
coordinates with dr = dO = 0, 


ds 2 = gudt 2 + 2 g t j, dt d(j) + 
or, after dividing by dt 2 , 


ds 


d(j) 


— ) -g tt + 2 g t(t> ^p+g H ( 


dt 


d(j) 


dt 


The explicit form of the metric components is not needed here. Note that 
ft = g is the angular velocity with respect to a distant observer, 

— 9tt + 2 gtcj, ft + g((,((, ft" . (143) 



The worldline of a particle has to be timelike, ds 2 < 0. Since the last equation 
is quadratic in ft, this is only possible between the roots ds 2 = 0, 


ft 


min / max 


9t(j> 

9<f>4> 



9tt 

9m 


What does 


f2 m i n ^ 


max 


mean?! In flat Minkowski spacetime (with Cartesian coordinates), ft m i n / max = 
±1 implies that a particle, e.g., may freely circle around a point, restricted 
only by the condition ]u| = |r • ft| < c. In the Kerr spacetime, at r = r E± , 
the smallest possible value of ft becomes 0. The particle may just stay at 
rest, but can rotate only in one direction, namely in direction of the angular 


56 



momentum of the black hole. Beyond te + , is forced to be larger than zero: 
The particle must co-rotate with the black hole. 

The preceding statement is only correct for radially infalling particles. In 
general, the influence of the rotating black hole on the motion of particles is 
more complex[165]. 

Rotation of the event horizon 

The behavior of on the event horizon is quite remarkable. By using the 
identities 

2 mr± = r± + a 2 , p 2 — r 2 = a 2 cos 2 6 = a 2 — a 2 sin 2 6 , (144) 

one finds for the event horizon (r = r + ), 

fl_|_ = := -. (145) 

2 mr + 

To interpret this result we use (143) and write 

Qtxv l* 1 l v \r=r+ — 0 i with P = (1, 0, 0, Q,h) . (146) 

The integral lines of l M = iT, 

aT = (t,r + ,9 0 ,tl H t), (147) 

define a lightlike hypersurface rotating with a uniform angular velocity: The 
event horizon of a Kerr black hole rotates “rigidly” with see in this 
context Frolov & Frolov[63]. A consequence of this finding is discussed in the 
next paragraph. 

Penrose process and black hole thermodynamics 

The (outer) ergosurface is a Killing horizon, not an event horizon. It is 
possible for particles to pass from the inside to the outside. This allows for 
a peculiar scenario: Since inside the Killing horizon the particle is forced to 
spin around, it picks up an additional rotational energy. This energy can 
be partly extracted by means of the Penrose process. An infalling particle 
traverses the Killing horizon, picks up rotational energy and subsequently 
decays into two parts. If one part plunges into the event horizon, the other 
part, carrying away some of the rotational energy, can return to the outside 
of the Killing horizon. Thus, the region between Killing and event horizon 
is justly labeled as “ergoregion” (from Greek ergon = work). 
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The observation that energy can also be extracted from the black hole 
gave rise to black hole thermodynamics. The next question is then how 
the parameters change if the black hole is infinitesimally disturbed. It was 
Bekenstein[13] who established a relation between the variations of the mass, 
the angular momentum, and the area of the event horizon. Using the coframe 
(185), for A = 0, we find for the area of the event horizon 



r = r_ i_ 
t = const. 


sin# (r 2 + a 2 ) ddAdcj) 


r = r_ i_ 
t = const. 


47t (r 2 + a 2 ) . (148) 


We can rewrite (148), using (144) and J = ma, 

A = 87rmr + = 87 t (m 2 — Vm 4 — J 2 ). 

The differential of this equation is 

dA d A 8n 87 t 

dA = —— dm + —— dJ = — dm - Uh dJ , 

dm OJKK 


with 


k = 


1 \Jm 4 — J 2 


2m rri 


+ y/m~‘ 


0/7 — K 


J 


(149) 


(150) 


(151) 


i4 — J 2 vm 4 — J 2 

The parameter Qjj is the angular velocity of the horizon (145). The param¬ 
eter k is the surface gravity. Eq.(150) can be rewritten as 


dm = — dA + O h dJ . (152) 

8tt 

The infinitesimal change of the mass, dm , is proportional to the the infinites¬ 
imal change of the energy, dE. The term O# dJ describes the infinitesimal 
change of the rotational energy. This suggests the identification of (152) 
with the first law of thermodynamics. The analogy is still more compelling 
by observing that, for a given black hole of initial (or irreducible) mass m, 
the area of the horizon is always increasing. Even by exercising a Penrose 
process, which extracts rotational energy from the black hole, a fragment of 
the incoming particle will fall into the black hole thereby increasing its mass 
and, in turn, the area of the horizon. Accordingly, the area A of the hori¬ 
zon behaves formally as if it is proportional to an entropy S and the surface 
gravity k as if it is proportional to a temperature T. In fact, the Hawking 
temperature and the Bekenstein-Hawking entropy turn out to be 


T 


h 

2nk-B 


«, 


S = 


AGh 


41, 


(153) 
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with k~B as the Boltzmann constant. Eq.(152) together with its thermody¬ 
namical interpretation (153) can be considerably generalized thereby estab¬ 
lishing the new discipline of “black hole thermodynamics”, see Heusler[85] 
and Car lip. [28] 


3.6 Beyond the horizons 

In the Schwarzschild spacetime, event horizon and Killing horizon coincide. 
In the Kerr spacetime, for m > a, there is an outer Killing horizon, an 
event horizon, an inner Killing horizon and an inner horizon. So far, all 
the coordinate systems we used for the Kerr metric show singularities at the 
outer and inner horizons r — r±. The construction of a regular coordinate 
system is possible along the same lines as for the Schwarzschild metric. Of 
course, the corresponding calculations are much more involved for the Kerr 
case. Therefore, we will give more a kind of heuristic approach to motivate 
Kruskal-like coordinates for the Kerr metric. 


Using light rays as coordinate lines 

Onr first task is to construct Eddington-Finkelstcin like coordinates for the 
Kerr metric by considering radial light rays. We restrict ourselves to the 
case 6 = 0 = (f. The Kerr metric in Boyer-Lindquist coordinates reduces to 
(0 = 0 — y p 2 = r 2 + a 2 = A + 2 mr): 

ds 2 = — ^ dt 2 + ^r- dr 2 . 
p 2 A 

Hence, for in-/out-going light rays, ds 2 = 0, we find 

o 2 r 2 i a 2 

dt = ± dr = ± --—-- dr 

A (r — r + )(r — r_) 


or, explicitly, 


it 



r 2 + a 2 r\ + a 2 r 2 _ + a 2 

——-- = r -1-In r — r+\ -In r 

r_)(r — r+) r + — r_ r + — r_ 


rO+const. 


(154) 

Unlike in the Schwarzschild spacetime, there form two event horizons, at 
r = r_ and r = r +1 respectively. However, as a —> 0, r_ goes to 0, whereas 
r .|_ approaches 2m and the Schwarzschild situation is reproduced. 

We next focus on the (Boyer-Lindquist) coordinates (t, r) and how the 
horizons etc. will appear in terms of the new coordinates. The other co¬ 
ordinates and the regularity of the metric is not addressed. However, all 
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the details can be found in the literature, see Refs. [20, 30, 78] . Using (154) 
analogously to (42), we introduce Eddington-Finkelstein like coordinates for 
Kerr, 


v := t + r + cr + In \r — r + \ — cm In |r — r_ 
u := t — r — cr_|_ In \r — r+| T cr_ In \r — r_ 

where (according to the notation in Ref. [20]) 

r± T a 2 mr± 

cr± := - = — . = . 

r+ - r_ V?n 2 - a 2 


(155) 

(156) 


(157) 


Again, we can get rid of the coordinate singularity by rescaling u and v 
analogously to (48). Since we have two horizons, r = r + and r = r_, we have 
to decide with respect to which singularity we rescale. We firstly choose r + 
and define, see (48), 


with 


v 


u 


•*(&)■ 

-»(-i 


r — r_i_ 2 r+» 

-e 2<t + 

r — r_ 2 


r — r+| 2 


r — r_ 2 


r—t 

e 2a + , 


(158) 

(159) 


i/:= — = — >1. (160) 

cr + T + 

Again, we go back to time- and space-like coordinates, exactly like in (50), 


t 


2 + “) > 


r : = 


1 

2 


(u — w) . 


(161) 


Then we work out the four coordinate patches exactly like (55) to (60). 
We arrive at a Kruskal like coordinate system, ffowever, there arises an 
important difference: The coordinate system still is singular for r = r_. This 
can be most easily seen from the analog to (59), the inverse transformation 
to r, which now reads 


~9 -9 ~ ~ r — r. 

r—t — —vu =-r— e ff + . 

(r - r_) v 


(162) 


The horizon r = r + is regular in this coordinate system and is described by 
r = ±t. The transformation(s) are valid in the domain r_ < r < +oo 


' + 


r = r_| 
r —> Too 
r —» r_ 
r = r E , 


r = ±t 

r 2 — t 2 —)• Too 


t 2 ^ 


-oo 


as for Schwarzschild 
particularly r —y Too for t — 0 
particularly t —y Too for r = 0 


t 2 = const. > 0 hyperbolas in I, 11 patches 
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In contrast to the Schwarzschild case, the full upper and lower halfplanes 
of the (r, t) plane is covered. It is not limited by the hyperbolas of the 
Schwarzschild singularity r = 0! 

We can regularize with respect to r_ by introducing 


v := — exp 




(163) 


Now we find 


P = 


T — r 

e CT - 


r+)' 


(164) 


This coordinate system covers the domain — oo < r < r + . Like the first 


coordinate system, it contains also the 
r < r + . This time, r >r + is excluded. 


r = 

: r 

'-to 

HH 

ii 

r — >- — oo 

: r 2 

— t ‘ 2 — > oo 

r —> r + 

: r 2 

1 

42 

i 

1 

8 

r = r E _ 

. r ,2 

— i 2 = const 

r = 0 

. f .2 

-t 2 = 


region between the horizons, r_ < 
as above 

particularly t —> ±oo for f = 0 
particularly r —>■ ±oo for t = 0 
hyperbola in I*, II* patches 
hyperbola in I*, II* patches 


Again, the whole (r, t ) plane is covered. Note that the spacetime extends 
beyond the ring(!) singularity. 


3.7 Penrose-Carter diagram and Cauchy horizon 

We compactify the Kruskal-like coordinate system for Kerr, yielding confor¬ 
mal Penrose-Carter diagrams. We discuss the analytical extension and the 
role of the inner horizon as Cauchy horizon. 

In order to draw a Penrose-Carter diagram for the Kerr spacetime, we 
compactify the coordinates via the tangent function like in Sec.2.6. The 
result looks at first quite similar to Schwarzschild in Fig. 11. However, the 
cutoff at r = 0 vanishes. The diagrams Fig. 13 and Fig. 14 both show the 
entire compactihed (r, f)-space. 

The two coordinate sets overlap in the region between the horizons. Thus, 
the corresponding coordinate patches have to be identified. And we can even 
draw beyond that ... Patch II is identified with patch IV*, II* with another 
patch IV**. And so on: We find an infinite sequence of coordinate systems. 
Formally, this constitutes a maximal analytic extension of the Kerr spacetime. 
Alas, there are good reasons for not believing in such vast an extension. 

The Kerr metric is a vacuum solution of Einstein’s held equation—it 
describes a totally empty spacetime. To render it physically meaningful, we 
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Figure 13: The Penrose diagram for the Kerr spacetime for r > r_. 
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Figure 14: The Penrose diagram for the Kerr spacetime for r < r_. 
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Figure 15: Maximal analytic extension of the Kerr spacetime 
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should regard it as the spacetime structure generated by a sensible physical 
source. One may ask then, why a single source should produce an infinite 
number of spacetimes. And it is even worse. The regions beyond the Cauchy 
horizon are exceptionally badly behaved. Consider the Cauchy surface in 
regions 1+11 of Fig. 15. All light rays and particle trajectories from the past 
intersect this surface only once. Then the held equations will tell us their 
future development, see Franzen[ !], for example. In Fig. 15 this is roughly 
indicated by the little light cone. However, even total knowledge of the world 
in I+II does not determine what might be going on in regions I*+II*. That is 
why r = r_ is called a Cauchy horizon, see Fig. 16. Thus, I* and II* are not 
only beyond the Cauchy horizon but also beyond predictable, sound physics. 
Moreover, the zigzagged region beyond the singularity is physically doubtful. 



Figure 16: Cauchy horizon: The causal past of a point P outside the Cauchy 
horizon of S is entirely determined by the information given on the Cauchy 
surface S. A point Q inside the Cauchy horizon receives also information from 
Z~. Evidently, initial data on S are not sufficient to uniquely determine events 
at point Q. The surface that separates the two regions “causally determined 
by S” and “not causally determined by S” is called Cauchy horizon. 

In this region, the asymptotics is reversed, see the permutation of I + and 
I~. As a consequence, the asymptotic mass in I* picks up a minus sign as 
compared to I. So the same source possesses a positive mass +m in I and a 
negative mass — m in I*, which seems strange. Moreover, it turns out that 
these regions are crowded with closed timelike curves. The whole extension 
is not globally hyperbolic. Thus one should restrict to the “diamond of sound 
physics”, I+II+III+IV. To do this consistently, one has to devise a physical 
mechanism preventing traveling beyond the Cauchy horizon, that is, the 
Cauchy horizon should become singular in some sense (cosmic censorship, 
see Penrose [150]). 
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3.8 Gravitoelectromagnetism, multipole moments 

The curvature tensor of the Kerr metric is calculated. By squaring it suitably, 
we find the two quadratic curvature invariants. Subsequently, we determine 
the gravitoelectric and the gravitomagnetic multipole moments of the Kerr 
metric, and we mention the Simon-Mars tensor the vanishing of which leads 
to the Kerr metric. 

The analogy between gravity and electrostatics became apparent when 
the Coulomb law was discovered in 1785. The gravitational and the elec¬ 
trostatic forces both obeyed an inverse-square law, with the difference that 
the mass can only be positive whereas the electric charge exists with both 
signs. Equal electric charges repel, opposite ones attract; in contrast, gravity 
is always attractive. 

In 1820 electromagnetism was discovered by Oersted, and the emerging 
unified theory, called “electrodynamics” by Ampere, eventually found its 
expression in the Maxwell equations of 1864. Besides the electric field E 
related to charge, we have the magnetic field B related to moving charge. 
These fields, together with the electric and magnetic excitations D and H , 
respectively, obey the Maxwell equations. 

Newton’s gravitational theory was only superseded in 1915/16 by Ein¬ 
stein’s gravitational theory, general relativity. However, already in the 1870s 
physicists began to speculate whether, besides Newton’s “gravitoelectric” 
field, related to mass at rest, there may also exist a new “gravitomag¬ 
netic” field, accompanying moving mass; for more details and references see 
Mashhoon[ .16]. As we saw above, these speculations became a solid basis 
in general relativity. In (90), the gravitomagnetic Lense-Thirring term sur¬ 
faced, which found solid experimental verification in the meantime. Thus, we 
can speak with justification of gravitoelectromagnetism[ ] (GEM), a notion 

which can guide our intuition, see in this context also Ni and Zimmermann[ 

Gravitoelectromagnetic field strength 

Electrodynamics is a linear theory, GR a nonlinear one. Still, if we take 
a linearized version of GR, there are those strong analogies between elec¬ 
trodynamics and gravitodynamics, as worked out, for instance, nicely in 
Rindler’s[165] book. However, the analogies go even further, as pointed out 
particularly by Mashhoon[ 16]. Even in an arbitrary gravitational field, if re¬ 
ferred to a Fermi propagated reference frame with coordinates (T, X), GEM 
is a useful concept. If we apply the geodesic deviation equation (33) to such a 
frame, the gravitoelectromagnetic field strength, representing the tidal forces, 
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turns out to be[116] 20 

G * M F a p = -R tm (T)X i . (165) 

If we develop (33) up to the order linear in the velocity V := dX/dT, we 
find 


d 2 X i 
dT 2 


-Roi ojX j + 2 R lk0j X j V k 


GEMp _ 2 GEMp y k 


(166) 


Now we recall that in electrodynamics the electric and the magnetic fields 
E and B, respectively, are accommodated in the 4d electromagnetic field 
strength tensor according to 


(F af ) 
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-E 1 

— E 2 

—E 3 
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-b 2 
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B 1 


V <> 
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= ~ (f>«) . 


(167) 


The diamond symbol o denotes matrix elements already known because of 
the antisymmetry of the matrix involved. The corresponding 2-form reads 
F = \F a p dx a A dx 13 . Keeping (167) in mind, Eq.(166) can be rewritten as a 
vector equation 

d 2 X 

— =-S r E-2Vx gr B. (168) 

dT 2 

In accordance with the equivalence principle, this equation of motion is in¬ 
dependent of the mass. The analogy with electromagnetism requires that the 
gravito electric charge, in terms of the mass m, is —1 and the gravito magnetic 
charge —2. In electrodynamics, both quantities are +1. The difference comes 
from the vector nature of the electromagnetic potential A a as compared to 
the tensor nature of the gravitational potential g a y, that is, hclicity 1 as 
compared to hclicity 2. The relation between the gravitomagnetic to the 
gravitoelectic charge, that is, the gyrogravitomagnetic ratio , is two: gr 7 = 2. 
Note that in Gravity Probe-B the authors specify the gyrogravitomagnetic 
ratio as 1. However, their gyros carried only orbital angular momentum 
rather than spin angular momentum. Hence this is to be expected; for more 
detailed discussions on this difference, see Ref. [80, 138]. 

20 Alternatively, we could generalize the Newtonian tidal force matrix of (9) to the grav- 
itoelectric and gravitomagnetic tidal force matrices, = Riojo and = CikiRkijo, 
respectively, see Scheel & Thorne[169]. Both matrices are symmetric and trace-free. Note 
that GEM F a p is an antisymmetric 4x4 matrix and £ and B are both symmetric trace-free 
3x3 matrices. 
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It has been pointed out by Ni[139] that the “measurement of the gyrograv- 
itational ratio of [a] particle would be a further step [138] towards probing 
the microscopic origin of gravity. GP-B serves as a starting point for the 
measurement of the gyrogravitational factor of particles.” 

Quadratic invariants 

In electrodynamics, we have two quadratic invariants [81]: 

^F af} F a P = *(*F A F) = B 2 - E 2 , ^F a pF* af3 = *(F A F) = 2E B , 

(169) 

where we used for the tensor dual the notation F* a P ■= i £ a/37 s F 7 g. We also 
employed the very concise notation of exterior calculus with the Hodge star 
operator. 21 

The first invariant is proportional to the Maxwell vacuum Lagrangian 
and is an ordinary scalar, whereas the second one corresponds to a surface 
term and is a pseudoscalar (negative parity). 

Turn now directly to the Kerr metric and list for this example the tidal 
gravitational forces, which are represented by the curvature tensor. With its 
20 independent components, it can be represented by a trace-free symmetric 
6x6 matrix, see (32). The collective indices A, B,.. = 1,..., 6 are defined as 
follows: {tf, t9, tif>\ 9(j) : (f>r, rQ} —> {1, 2,3; 4, 5, 6}. We throw the orthonor¬ 
mal Kerr coframe (128) to (133) into our computer and out pops the 6x6 
curvature matrix, 


with 
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ma cos 

e 3r2 


= ( Rba ) , 


3r 2 — a 2 cos 2 6 
r 2 + a 2 cos 2 6 ) 3 


(170) 


(171) 


It is straightforward to identify E as the gravitoelectric and B as the gravit- 
omagnetic component of the curvature. This is in accordance with (165). 

21 The Hodge star of a p-form w = (l/p\)u} fJ , 1 ... f _ lp dx IJ ’ 1 A • • - A dx^ p is an [n— p)-form *w, 
with the components (*w) Ml ... M „_ p = (l/p!)e l/1 where e is the totally 

antisymmetric unit tensor and n the dimension of the space, see Eq.(C.2.90) in Ref. [8 ] 
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It is obvious how we should continue. Our gravitoelectromagnetic invari¬ 
ants will be 22 

K := ^R a ^ s R a0lS = ~ *{*R«p A R a ?) , (172) 

V := = . (173) 

Again, our program determines the Kretschmann 23 scalar K and the Chern- 
Pontryagin pseudoscalar V to be[ .9] 

K = —24(1B 2 - E 2 ), V = —48 EB . (174) 

The similarity to (169) is impressive. The GEM analogy quite apparently 
applies to the full nonlinear theory. The results in (174), partly in more 
involved representations, can be found in the literature, see, for instance, the 
books of de Felice & Clarke[42] and of Ciufolini & Wheeler[36], but compare 
also de Felice & Bradley) II], Henry) 1], and Cherubini et al.[ ]. 

Thus, the quadratic invariants K and V confirm that the Kerr metric 
is the exterior field of a rotating mass distribution. In order to get more 
information about this distribution, we proceed, like in electrodynamics, and 
look into the gravitoelectromagnetic multipole moments of this rotating mass. 

Gravitomagnetic clock effect of Mashhoon, Cohen, et al. 

According to the results of Lense-Thirring, the rotation of the Sun changes 
the spacetime around it by inducing gravitomagnetic effects. As we saw 
above, in a similar way the temporal structure around a Kerr metric is af¬ 
fected by the angular momentum of the Kerr source. Thus, a gravitomagnetic 
clock effect should emerge, 24 the measurability of which requires very accu¬ 
rate clocks. The effect can be demonstrated by two clocks that move on 
equatorial orbits, one in prograde and the other in retrograde orbit around 
the Kerr metric. It turns out [117] that the prograde equatorial clock is slower 
than the retrograde one. This is not necessarily what our intuition would 

22 In exterior calculus, we have the Euler 4-form E := R a p A with K = *E. 

Analogously, we have the Chern-Pontryagin 4-form P := —R a p A R^ a , which is an exact 
form, with V := *P, cf. Obukhov et al. [141]. 

23 Usually in the literature [36, 42], the Kretschmann scalar is defined as R a / 3 ^ 5 R al3lS , 
even though the electrodynamics analogy would suggest to include the factor 1/2. 

24 This was first predicted by Cohen and Mashhoon[37] and worked out in greater detail 
by Mashhoon et al. [118, 117], see also Bonnor & Steadman[18] and the review papers in the 
workshop of Lammerzahl et al. [103] . In a similar way, there emerges also a gravitomagnetic 
time delay, see Ciufolini et al. [33] . 
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tell us. It is connected with the fact that the dragging of frames in a Kerr 
metric can sometimes turn out to be an “antidragging”, thus making this 
notion less intuitive, [165] as we already recognized in Sec.3.5. 

Generalizations of this clock effect were studied, for example, by Hack- 
mann & Lammerzahl[ ]. The recent discussion of the Clocks around SgrA*, 
by Angclil & Saha[l] is, in effect, just one more manifestation of the gravit- 
omagnetic clock effect. 

Multipole moments: gravitoelectric and gravitomagnetic ones 

In Newton’s theory, one gets a good idea about a mass distribution and 
its gravitational field by determining the multipole moments of the mass 
distribution M. In GR, because of the existence gravitomagnetism, we have 
to expect a new type of multipole moments, namely the moments J of the 
angular momentum distribution. 

If a stationary axially symmetric line element of the form (93) is asymp¬ 
totically flat, then it is possible[L80] to define two sets of multipole mo¬ 
ments, the gravitoelectric moments M s (“mass multipole moments”) and the 
gravitomagnetic moments J s (“angular momentum multipole moments”), for 
s = 0,1, 2,.... These moments were found by Geroch[ ] for the static and by 
Hansen[ ] for the stationary case. They were reviewed by Quevedo[160] and 
used for constructing new exact solutions by Quevedo & Mashhoon[ L59, 161]. 
Hansen computed the multipole moments for the Kerr solution and found 


s = 0 

M 0 = 

—m 

Ji = ma 

(175) 

s = 1 

m 2 = 

ma 2 

J 3 = —ma 3 

(176) 

s = 2 

M 4 = 

4 

—ma 

J .5 = ma 5 

(177) 

s = 3 ... 

m 6 = 

ma 6 .. 

= —ma 7 ... 

(178) 

More compactly, 

we have 




m 2s = (- 

l) s+1 ma 2s , 


M 2s +i = 0; 

(179) 

— 0 , 



J 2s +i = (—1 ) s ma 2s+1 . 

(ISO) 


It is possible to introduce normalized multipole moments, see Meinel et 
al.[ 3], such that for Kerr we have M s + iJ s = m(ia) s . Then the mass 
monopole Mo = m is positive. Apparently, the Kerr metric has a simple 
multipolar structure or, formulated differently, only very specific matter dis¬ 
tributions can represent the interior of the Kerr metric. 

Quevedo[L60] compiled a number of theorems which illustrate the use of 
the multipole moments: 
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1) A stationary spacetime is static if and only if all its gravito magnetic 
multipole moments vanish (Xanthopoulos 1979). 

2) A static metric is flat if and only if all its gravito electric multipole 
moments vanish (Xanthopoulos 1979). 

3) A stationary metric is axisymmetric if and only if all its multipole 
moments are axisymmetric (Giirsel 1983). 

4) Two metrics with the same multipole moments have the same geometry 
at large distances from the source (Beig & Simon 1981; Kundu 1981; 
Van den Bergh & Wils 1985). 

5) Any stationary, axisymmetric, asymptotically flat solution of Einstein’s 
vacuum equation approaches the Kerr solution asymptotically (Beig & 
Simon 1980). 

6) Any static, axisymmetric, asymptotically flat vacuum solution approaches 
the Schwarzschild solution asymptotically (Beig 1980). 

In the formulation of Stephani, Kramer, et al. [180]: 

7) A given asymptotically flat stationary vacuum spacetime is uniquely 
characterized by its multiple moments. 

We recognize that the knowledge of the multipole moments provides a lot of 
insight into the physical properties of an exact solution. 

From the point of view of the Kerr solution, theorem 5), see Beig & 
Simon[12], is perhaps the most interesting one. It underlines the central 
importance of the Kerr solution. The considerations in the context of theorem 
5) were further developed by Simon[ 7 4, 17 ]. On the 3-dimensional spatial 
slices of a stationary axially symmetric metric, he defined the 3d “Simon 
tensor”,[ 15] a kind of complexified generalized Cotton-Bach tensor[65]. The 
vanishing of the Simon tensor then leads to the multipole moments of the Kerr 
solution. Later, Mars[112], see also Mars[: ] and Mars & Senovilla[113], 
generalized this approach and was led to the 4d “Simon-Mars tensor”. In 
Ionescu & Klainerman[88], one can find a more extended discussion of the 
Simon-Mars tensor, see also Wong[192], More recently, Backdahl & Valiente 
Kroon[9] have proposed replacing the Simon-Mars tensor by another measure 
of “non-Kerrness”, namely a scalar parameter. 
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3.9 Adding electric charge and the cosmological con¬ 
stant: Kerr-Newman metric 

Enriching the Kerr metric by an electric charge is straightforwardly possible. 
We start from the metric (133) with coframe (128) to (132). This coframe 
can accommodate the Kerr, the Schwarzschild, and the Reissner-Nordstrom 
solutions. The different forms of the function A suggest how a charged Kerr 
solution should look like ... 


Schwarzschild 

(m) 

P = 

r 2 

A = 

r 2 

— 2 mr 


Reissner-Nord. 

{m, q ) 

P = 

r 2 

A = 

r 2 

— 2 mr 

+q ‘ 2 

Kerr 

(m, a) 

P = 

r 2 + a 2 cos 2 6 

A = 

r 2 

— 2 mr 

+a 2 

Kerr-Newman 

(m, a, q) 

P = 

r 2 + a 2 cos 2 9 

A = 

r 2 

— 2 mr 

+q 2 +a 2 


Charging the Schwarzschild solution is achieved by adding q 2 to the func¬ 
tion A. Since the charged Kerr solution should encompass the Reissner- 
Nordstrom solution, we tentatively keep the term q 2 for the case a ^ 0. 
Now, we can indeed find a potential, 


A 



(dt — a sin 2 6 d(j >), 


(181) 


such that the Einstein-Maxwell equations are fulfilled. The potential de¬ 
scribes a line-like charge distribution at p — 0, that is, on the ring singu¬ 
larity of the Kerr spacetime, which is quite satisfying[135] . This charged 
Kerr solution was first worked out by Newman, Couch, Chinnapared, Ex¬ 
ton, Prakash, and R. Torrence[134] (1965), using “methods which transcend 
logic”, as Ernst[54] puts it. He, in turn, proceeded from (120). Replacing 25 
£ by yR — qq* £ generates a solution of the Einstein-Maxwell equations with 
potential A t + iA ^ = <?/(£ + 1). 

The Kerr and the Kerr-Newman solution behave quite similarly. We can 
adopt most of the discussion of the Kerr metric by substituting a 2 + q 2 for 
a 2 . 

We can further generalize the Kerr-Newman metric to include also a cos¬ 
mological constant, see Sec.4.1., and even more parameters, see Fig. 17. 

25 Here, q is not the charge but a complex parameter in the solution of the Ernst equation. 


71 






Plebanski —Demianski 



Figure 17: Schematics of Petrov D solutions: The spherically symmetric 
Schwarzschild solution with mass parameter m is located in the center. 
Adding an electric charge q brings us to the Reissner-Nordstrom solution. 
It is still spherically symmetric but adds a second horizon. The distance 
between the horizons increases with the charge q. Setting the black hole into 
rotation, the angular momentum parameter emerges, a ^ 0, and reduces the 
spherical symmetry to an axial one. An oblate ergosurface (two, actually) 
forms. Event horizon and ergosurface meet at the polar axis, the equatorial 
distance increases with a. All these solutions can be deSittere d, that is, a 
cosmological constant A is added. All presented solutions are subcases of the 
Plebanski-Demianski solution, which adds three more parameters [69]. 


3.10 On the uniqueness of the Kerr black hole 

The Kerr black hole, up to some technical assumptions, is the unique solution 
for the stationary, axially symmetric case. We point to some of the literature 
where these results can be found. 

Because of the Birkhoff theorem, the Schwarzschild solution (mass param¬ 
eter m ) represents the general spherically symmetric solution of the Einstein 
vacuum field equation. The analogous is true in the Einstein-Maxwell case 
for the 2 parameter Reissner-Nordstrom solution (mass and charge parame¬ 
ters m and q, respectively). Thus, for spherical symmetry, we have a fairly 
simple situation. 
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In contrast, in the axially symmetric case, there does not exist a gener¬ 
alized Birkhoff theorem. The 2-parameter Kerr solution (mass and rotation 
parameters m and a, respectively), is just a particular solution for the axi¬ 
ally symmetric case. As we saw in Sec.3.8, the Kerr solution has very simple 
gravitoelectric and gravitomagnetic multipole moments (179,180). Numer¬ 
ous solutions are known that represent the exterior of matter distributions 
with different multipole moments. The analogous is valid for the 3 parame¬ 
ter Kerr-Newman solution (parameters m,a,q ), see Stephani et al. [180] and 
Griffiths & Podolsky [69]. 

However, one can show under quite general conditions that the Kerr- 
Newman metric represents the most general asymptotically flat, stationary 
electro-vacuum black hole solution (“no-hair theorem”), see Meinel’s short 
review[ ]. Important contributions to the subject of black hole uniqueness 
were originally made by Israel[90, 91], Carter[30, 31], Hawking & Ellis [78], 
Robinson[166, 167], and Mazur [120] (1967-1982), for details see the recent 
review of Chrusciel et al. [ ]. 

More recently Neugebauer and Meinel[130, 133] found a constructive 
method for proving the uniqueness theorem for the Kerr black hole met¬ 
ric. This was extended to the Kerr-Newman case by Meinel.[12 ] By inverse 
scattering techniques, they showed how one can construct the Ernst potential 
of the Kerr(-Newman) solution amongst the asymptotically flat, stationary, 
and axially symmetric (electro-)vacuum spacetimes surrounding a connected 
Killing horizon. 

Let us then eventually pose the following questions [2 7]: 

(i) Are axially symmetric, stationary vacuum solutions outside some mat¬ 
ter distribution “Kerr”? The answer is “certainly not”, and it makes 
sense to figure out ways to characterize the Kerr metric, see Sec.3.8. 

(ii) Is the Kerr solution the unique axially symmetric, stationary vacuum 
black hole? The answer is essentially “yes” (modulo some technical 
issues)—see, for example Mazur[120]. 

The general tendency in the recent development of the subject is to use 
additional scalar or other matter fields. They weaken the uniqueness theo¬ 
rems, which is probably not too surprising. 

Let us conclude with a quotation that may make you curious to learn 
still more about the beauty of the Kerr metric: We have many different 
axially symmetric solutions. The Kerr solution is characterized by “station¬ 
ary, axially symmetric, asymptotically flat, Petrov type D vacuum solution of 
the vanishing of the Simon tensor, admitting a rank-2 Killing-Stackel (KS) 
tensor of Segre type [(11)(11)J constructed from a (non-degenerate) rank-2 
Killing-Yano (KY) tensor ”, see Hinoui et al. [87] . 
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3.11 On interior solutions with material sources 

To match the Kerr (vacuum) metric to a material source consistently is one 
of the big unsolved problems. Only the rotating disc solution of Neugebauer 
& Meinel provides some hope. 

This section is added in order to draw your attention to an unsolved prob¬ 
lem, to the solution of which you might want to contribute. Find a realistic 
material source for the Kerr metric in the sense of an exact solution. Many 
unsuccessful attempts have been made, see the early review of Krasinski[101] 
of 1978. More recently, in 2006, Krasinski[ 56] concludes “that a bright new 
idea is needed, as opposed to routine standard tricks tested so far. ” This 
statement was not made lightheartedly, Krasinski knows what he is talking 
about. 

Many axially symmetric vacuum solutions were constructed. Quevedo & 
Mashhoonfl .61], for example, deformed the multipole moments of the Kerr(- 
Newman) metric and constructed appropriate solutions of the Einstein(- 
Maxwell) equation that describe the exterior gravitational field of a (charged) 
rotating mass. It is always the hope that somebody may find a suitable mat¬ 
ter distribution with the multipole moments of the Kerr solution—but this 
did not happen so far; for another approach see Marsh [ ]. 

We are only aware of one exact solution that fits into this general context: 
It is the infinitesimally thin and rigidly rotating dust solution of Neugebauer 
& Meincl[131, 132] (1993). It is an exact analytical solution of the Einstein 
equation with matter. It depends on 2 independent parameters, the radius 
p 0 of the disk and its angular velocity 12. Petroff & Meinel[151] developed, 
by means of an iterative procedure, a post-Newtonian approximation of the 
solution that helps to understand the Newtonian limit. 

We recall that in electrostatics in flat space, for example, we prescribe 
an electric charge distribution and we are used to solve the corresponding 
boundary value problem within Maxwell’s theory. Similarly, Neugebauer & 
Meinel specified a very thin rotating disk of dust and solved the boundary 
value problem within GR. This is a well-defined procedure. The problem is, 
however, that within a non-linear theory, such as GR, it is extremely hard 
to implement. Remarkably, for certain parameter values, the gravitational 
field of the disk approach the extremal Kerr case. Accordingly, there exists 
a certain relation to the Kerr problem. The desideratum would be to find 
a rotating matter distribution the external field of which coincides with the 
complete Kerr field. 

Driven by the fact that the electrically charged Kerr solution, the Kerr- 
Newman solution, has a g-factor of 2, exactly like the electron (see also Pfister 
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& King[153]), Burinskii[24, 25] speculated that a soliton like solution of the 
Dirac equation may be the source of the Kerr metric, see also Burinskii & 
Kerr [26]. Is that the “bright new idea” Krasinski was talking about? We do 
not know but a hard check of the Burinskii ansatz seems worthwhile. 


4 Kerr beyond Einstein 

In generalizations of Einstein’s theory of gravity, the Riemannian geometry 
of spacetime is often extended to a more general geometrical framework. We 
describe two such examples in which the Kerr metric still plays a vital role. 

4.1 Kerr metric accompanied by a propagating linear 
connection 

We display the Kerr metric with cosmological constant that, together with 
an explicitly specified torsion, represents an exact vacuum solution of the 
two field equations of the Poincare gauge theory of gravity with quadratic 
Lagrangian. 

In gauge theories of gravitation, see Blagojevic et al. [ 7], the linear con¬ 
nection becomes a field that is at least partially independent from the met¬ 
ric. It can be either metric-compatible, then it is a connection with values in 
the Lie-algebra of the Lorentz group 5'0(1,3) and the geometry is called a 
Riemann-Cartan geometry, or it can be totally independent, then it resides 
in a so-called metric-affine space and the connection is GL( 4, i?)-valued. For 
simplicity, we concentrate here on the former case, the Poincare gauge theory 
of gravity, but the latter case is also treated in the literature [ S7, 11]. 

Let us shortly sketch the theory. Gauging the Poincare group leads to a 
spacetime with torsion T a and curvature R a13 (Riemann-Cartan geometry 26 ): 

T a := Dd a = dd a + Tg a A 0? = ^T ij a dx i A dx j , (182) 

R af) dr af3 _ r a 7 A r ^/3 = _ R fta = A d X j . (183) 

Besides the coframe 1-form d Q , the Lorentz connection 1-form T a/3 = Ti a/3 dx z = 

_ T p a 

is a second field variable of the gauge theory. For a Riemannian space, 
torsion T a — 0 and T al3 becomes the Levi-Civita connection. 

26 Experimental limits of a possible torsion of spacetime were recently specified in a 
remarkable paper by Obukhov et al., [142] see also the literature given there. 
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We choose a model Lagrangian quadratic in torsion and curvature, in 
actual fact (for h — 1, c = 1), 


V 


1 

2k 


(T a A #) A *(7> A d a ) 


1 

2 ^ 


R a ^ A -kR a p , 


(184) 


with Einstein’s gravitational constant k (dimension length-squared) and a 
dimensionless strong gravity coupling constant g. One can calculate the two 
vacuum field equations by varying with respect to d a and T a/3 . In 1988, for 
these two field equations, a Kerr metric with torsion[10] was found as an 
exact solution. 

We display here the orthonormal coframe and the torsion: The coframe 
in terms of Boyer-Lindquist coordinates (t,r,6,(j)), reads (in the conven¬ 
tions used in Ref. [10]) 


( dt + a sin 2 0 dcf) , 

(185) 


(186) 

Jr*- 

(187) 

\fF sin 9 r , , 9 9n , n 

- [adt + (r 2 + a") d(f>\ . 

(188) 


As before, we have p 2 := r 2 + a 2 cos 2 9. However, the other structure functions 
pick up a cosmological constant A: 


F 1 + -Aa 2 cos 2 9 , A := r 2 + a 2 — 2 Mr — -Ar 2 (r 2 + a 2 ). (189) 

3 3 

The corresponding metric is called a Kerr-deSitter metric. The coframe is 
orthonormal. Then the metric reads 


g = —{P 0 0 t?* + t? 2 0 d 2 + 0 iP . (190) 


It is a characteristic feature of these exact solutions that even though the 
Lagrangian (184) does not carry a cosmological constant, in the coframe and 
the metric there emerges such a constant, namely A := —3q/(4k). This could 
be of potential importance for cosmology. 
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The torsion T a of this stationary axially symmetric solution of the Poincare 
gauge theory reads ($ a/3 := d a A P' 9 ), 
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(191) 


with the following gravitoelectric and gravitomagnetic functions: 
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( 192 ) 

( 193 ) 


Metric and torsion of this exact solution are closely interwoven. Note, in 
particular, that the leading gravitoelectric part in the torsion, for small a, 
is ~ M/r 2 , a definitive Coulombic behavior proportional to the mass. For 
a = 0, we find a Schwarzschild-deSitter solution with torsion. 

One may legitimately ask, why is it that the Lagrangian (184) yields 
an exact solution with a Kerr-deSitter metric? The answer is simple: The 
Lagrangian was devised such that the torsion square-piece, in lowest order 
in ft, encompasses a Newtonian approximation. This is already sufficient in 
order to enable the existence of a Kerr-deSitter metric. One could even add 
another torsion-square piece to V for getting an Einsteinian approximation, 
but this is not even necessary. Thus, only a Newtonian limit of some kind 
seems necessary for the emergence of the Kerr structure. 


4.2 Kerr metric in higher dimensions and in string the¬ 
ory 

There exist also Schwarzschild and Kerr metrics in higher dimensional space- 
times. These investigations are mainly motivated by supergravity and string 
theory. 

Tangherlini[ ] (1963) started to investigate higher dimensional Schwarzschild 
solutions, with n — 1 spatial dimensions. He studied the (“planetary”) or¬ 
bits in an n-dimensional Schwarzschild field (“Sun”) and found that only for 
n = 4 we have stable orbits, see also Ortin[ ]. According to Tangherlini, 
this is then the only case that is interesting for physics. Nowadays, however, 
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many physicists hypothesize that higher dimensions do exist because string 
theory suggests it. 

Somewhat later, Myers and Perry [128] (1986) generalized these consid¬ 
erations to higher-dimensional Kerr metrics. In the meantime a plethora 
of such higher-dimensional objects have been found, see Allahverdizadeh et 
al. [ ] and Frolov & Zelnikov[ ]. Recently Keeler et al. [97] investigated, in 
the context of string theory, the separability of Klein-Gordon or Dirac fields 
on top of a higher-dimensional Kerr type solutions. Lately Brihaye et al. [ 1], 
for example, discussed the exact solution of a 5d Myers-Perry black holes as 
coupled to a to a massive scalar held. The physical interpretations of these 
results remain to be seen. 
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Appendix 

A Exterior calculus and computer algebra 

We want to use as input the Papapetrou metric (93). We take the equivalent 
representation in the form of the orthonormal coframe of the Eqs.(94) to 
(98). How such a Reduce-Excalc program can be set up, is demonstrated 
in Stauffer et al. [178] and in Socorro et al. [ 76], for the Einstein 3-form, see 
Heinicke[82]: 

% Coframe of Andress-Lewis-Papapetrou-Buchdahl metric 

% file Buchdahl03.exi, 29 July 2014, fwh & chh 
% in "Buchdahl03.exi"; 

load_package excalc; 
off exp$ 

pform f=0, omega=0, garama=0 $ 

fdoraain f=f(rho,z), omega=omega(rho,z), gamma=gamraa(rho,z) ; 

coframe o(0) = sqrt(f) * (d t - omega * d phi), 

o(l) = sqrt(f)**(-1) * exp(gamraa) * d rho, 
o(2) = sqrt(f)**(-1) * exp(gamma) * d z, 
o(3) = sqrt(f)**(-1) * rho * d phi 
with signature (1,-1,-1,-1); 

displayframe; 
frame e$ 

% Connection, curvature, and Einstein forms 

pform connl(a,b)=l, curv2(a,b)=2$ 
antisymmetric connl, curv2$ 
factor o(0), o(l), o(2), o(3)$ 

connl(-a,-b) := (1/2)*( e(-a)_|d o(-b) - e(-b)_|d o(-a) 

- (e(-a)_|(e(—b)_|d o(-c))) * o(c))$ 
curv2(-a,b) := d connl(-a,b) - connl(-a,c) ~ connl(-c,b)$ 
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% Einstein tensor = Einstein 0-form 
pform einstein3(a)=3, einsteinO(a,b)=0$ 
symmetric einsteinO$ 

einstein3(-a) := -(1/2) * curv2(b,-c) ~ # (o(-a) ~ o(-b) ~ o(c))$ 
einsteinO(a,-b):= #( o(a) ~ einstein3(-b))$ 

on exp, gcd$ 
factor ~$ 
on nero; 

einsteinO(a,-b):= #( o(a) ~ einstein3(-b)); 
off nero; 

% by inspection, we find 

einsteinO(l,-l) + einstein0(2,-2) ; °/„ equals 0 
einstein0(0,-0) - einstein0(3,-3) ; °/ 0 eliminates gamma 

out "Buchdahl03.exo"; 

load_package tri; 
on tex; 
on TeXBreak; 

einsteinO(a,-b):=einsteinO(a,-b); 
off tex; 

einsteinO(a,-b):=einsteinO(a,-b); 
omega:=0; 

einsteinO(a,-b):=einsteinO(a,-b); 

shut "Buchdahl03.exo"; 

; end; 
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